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Important progress in understanding the behavior of hadronic matter at high den- 
sity has been achieved recently, by adapting the techniques of condensed matter 
theory. At asymptotic densities, the combination of asymptotic freedom and BCS 
theory make a rigorous analysis possible. New phases of matter with remarkable 
properties are predicted. They provide a theoretical laboratory within which chiral 
symmetry breaking and confinement can be studied at weak coupling. They may 
also play a role in the description of neutron star interiors. We discuss the phase 
diagram of QCD as a function of temperature and density, and close with a look 
at possible astrophysical signatures. 
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1 Introduction and Summary 

In this article we shall discuss the behavior of QCD at high density. We shall 
mainly be concerned, more precisely, with the regime of very large baryon 
number density and relatively low temperature. The appeal of the subject 
should be obvious. It provides the answer to a child-like question: What 
happens to matter, as you squeeze it harder and harder? Moreover, this regime 
may be realized in neutron star interiors and in the violent events associated 
with collapse of massive stars or collisions of neutron stars, so it is important 
for astrophysics. Finally, we may hope to gain insight into QCD at moderate 
or low density — the realm of ordinary matter and of terrestrial experiments 
— by extrapolating from above. 

What is less obvious, but turns out (with certain qualifications) to be true, 
is that the high density regime of QCD is tractable Iheoretically. Heuristic 
arguments to this effect, invoking asymptotic freedomP go back to the earliest 
days of the subjectcl High density brings in a large energy scale, the chemical 
potential /i, and one might hope that the relevant coupling to describe the 
dynamics is g{pL). This becomes small as /i — *■ oo, and suggests the possibility 
of a perturbative treatment. This naive expectation, however, does not stand 
up to critical scrutiny. As we shall discuss at length below, perturbation theory 
around the naive ground state (free quark Fermi spheres) encounters infrared 
divergences. Furthermore, the naive perturbative ground state is unstable. 
Therefore, straightforward perturbative treatment of QCD at high density fails. 

Fortunately, related difhculties have been met and overcome previously, in 
the theory of superconductivity. There we learn that arbitrarily weak attrac- 
tive interactions can change the ground state qualitatively. In the true ground 
state there is an effective mass for photons — the Meissner effect — and en- 
ergy gaps for charged excitations. These phenomena remove potential infrared 
divergences, and render the perturbation theory around the true ground state 
regular (nondegenerate). 

We can readily adapt the methods of superconductivity theory to QCD. It 
is also instructive to consider, in this connection, variants of QCD with different 
numbers of colors and flavors, and different spectra of quark masses. A rich 
and highly structured theory emerges, displaying calculable, highly non-trivial 
dependence on all these variables. 

The central result of the analysis is the identification of condensates in 
diquark channels, analogous to the Cooper pairs of electrons in ordinary su- 
perconductors. This is the phenomenon of color superconductivitydU'Era 

Compared to ordinary superconductivity, color superconductivity, though 
it appears superficially to be more complex mathematically, is in a profound 
sense simpler and more directly related to fundamentals. Ordinary super- 
conductivity takes place in solids and the accurate effective interactions are 
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determined by band structure and other complicated effects. Furthermore, 
ordinary superconductivity in a metal involves electron pairing, and the fun- 
damental interaction between electrons (the screened Coulomb interaction) is 
repulsive. The effective attraction near the Fermi surface that leads to su- 
perconductivity arises in classic superconductors only as a subtle consequence 
of retarded phonon interactions, and in the cuprate superconductors through 
some mechanism yet unknown. In color superconductivity, by contrast, the 
attractive interaction can arise already from the primary, strong, interactions. 
This has two consequences. First, the accurate form of these interactions can 
be calculated from first principles, using asymptotic freedom. This makes 
calculations at high enough density robust. Second, at accessible densities, 
where the strong interactions are much stronger than the electromagnetic in- 
teractions, we expect the color superconductors themselves to be robust in the 
sense that the ratio of their gaps and critical temperatures to the Fermi energy 
should be quite large. 

In QCD with three colors and three flavors, we find an improved ground 
state at high density, based on color superconductivity, around which weak- 
coupling perturbation theory is valid. In particular, all the colored degrees of 
freedom acquire gaps. Thus, the improved ground state differs qualitatively 
from the naive one. 

The resulting predictions regarding the low-energy spectrum and dynamics 
are striking. Color symmetry and chiral symmetry are spontaneously broken. 
The spectrum of elementary excitations is quite different from that found in 
naive perturbation theory. Nominally massless quarks and gluons become mas- 
sive, new massless collective modes appear, and various quantum numbers get 
modified. All the elementary excitations carry integral electric chargesB Alto- 
gether, one finds an uncanny resemblance between the properties one computes 
at asymptotic densities, directly from the microscopic Lagrangian, and the 
properties one expects to hold at low density, based on the known phenomenol- 
ogy of hadrons. In particular, the traditional "mysteries" of confinement and 
chiral symmetry breaking are fully embodied in a controlled, fully microscopic, 
weak-coupling (but nonperturbative!) calculation—that accurately describes a 
directly physical, intrinsically interesting regimeOu 

Though some of the ideas go back a decade or more^Q the full power of 
color superconductivity to provide a rigorous foundation for the investigation 
of high density QCD has only become apparent relatively recently, and the 
subject is developing rapidly. In this survey we shall emphasize what we see 
as the most fundamental ideas that have appeared in the field to date, and 
attempt to identify some significant challenges for the future. Although we 
shall supply extensive references we will not attempt to catalogue all the very 
latest results, nor snapshot the developing state of the art in its technical and 
quantitative aspects. 
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In Section 2, we shall discuss in detail the high-density behavior of a 

slightly idealized version of real-world QCD, in which we imagine there are 
three flavors of massless quarks. This is the case for which the clearest and 
most beautiful results emerge, based on the phenomenon of color-flavor lock- 
ing. In Section 3, we shall discuss what changes must be made, in order to 
take into account the realistic quark spectrum. We shall also briefly discuss 
QCD variants with unrealistic numbers of colors and flavors, which shed addi- 
tional light on the theory. In Section 4, we shall discuss in greater depth the 
theoretical analysis which underlies the preceding sections, and the techniques 
it relies upon. In Section 5, we shall synthesize the preceding discussion, to 
produce a tentative sketch of the phase diagram for high-density QCD. Finally 
in Section 6 we shall discuss possible applications to astrophysics, including 
both existing results and some avenues that we believe might reward further 
investigation. 

2 Color-Flavor Locking 

In this Section we shall analyze the high-density, zero-temperature behavior 
of a slight idealization of QCD, in which the masses of the u, d and s quarks 
arc set to zero, and those of the c, b and t quarks to infinity. This idealization 
gives rise to an especially clear and beautiful form of the theory. Also, as we 
shall discuss in Sections 3 and 5, the analysis applies with only very minor 
modifications to an important phase of real- world QCD. 

In this Section our discussion will be broadly conceptual. The formal and 
algorithmic underpinnings are spelled out in more detail in Section 4, and of 
course in the primary literature. In particular, wc focus in this section on 
a presentation of the physical properties of dense quark matter in the ideal- 
ized three-flavor world, deferring discussion of quantitative calculations of the 
magnitude of the gap at the Fermi surface as much as possible to Section 4. 

Let us briefly describe the foundational argument justifying the weak- 
coupling approximation at high density, leaving a more detailed and precise 
discussion to Section 4. The relevant degrees of freedom in cold, dense quark 
matter are those which involve quarks with momenta near the Fermi surface. 
At high density, when the Fermi momentum is large, the QCD gauge coupling 
g{p) is small. When a pair of quasi-particles scatter, the typical momentum 
transfer is of order /x and the interaction is therefore weak. The exception, 
of course, is the case of scattering by a small angle. In QCD at asymptoti- 
cally high densities, the long-range magnetic gauge interactions are unscreened 
in the absence of superconductivity, which raises the possibility of infrared 
problems in small angle scattering. The Meissner effect induced by the su- 
perconducting condensate itself can provide the requisite screening, however, 
regulating the putative coUinear divergence and guaranteeing that the calcula- 
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tion is truly a weak-coupling calculation, at least at truly asymptotic density. 
Quantitative calculations support this qualitative conclusion, although it turns 
out that dynamical screening is more important than the Meissner effect in the 
regulation of small angle scattering. 

2. 1 Form of the Condensate 

Because of the infinite degeneracy among pairs of quarks with equal and oppo- 
site momenta at the Fermi surface, an attractive interaction between quarks, 
even if arbitrarily weak, renders the Fermi surface unstable to the formation 
of a condensate of quark Cooper pairs. Creating a pair costs no free energy at 
the Fermi surface, while the attractive interaction results in a free energy bene- 
fit. The consequence of an attractive interaction between opposite-momentum 
modes near the Fermi surface is therefore the formation of a condensate in the 
zero temperature ground state. We expect those quark quasiparticles which 
interact with the condensate to acquire an energy gap, and we expect a Meiss- 
ner effect to occur for all gauge bosons except those which see the condensate 
as neutral. 

Single gluon exchange, the QCD analogue of the Coulomb interaction be- 
tween two quarks, is attractive if the quarks are antisymmetric in color and 
the pair is therefore in the color 3 channel. At weak-coupling, this interaction 
dominates and this argument suffices to guarantee condensation in the color 3 
channel. The instanton interaction is also attractive in the 3 channel, which 
may be relevant at stronger coupling. At any coupling, attraction in the 3 
channel is quite reasonable intuitively, for in this channel the total color flux 
is reduced as one brings the quarks together. 

Thus, according to the preceding discussion, we should expect the forma- 
tion of a condensate in a color 3 channel. The form of the condensate can be 
anticipated by analyzing the renormalization of couplings toward the Fermi 
surface starting from realistic microscopic couplings, as we describe in Sec- 
tion 4. This indicates the most likely flavor and spin channels for condensate 
formation. Ultimately, one must compare the energies of different candidate 
ground states, constructed along the Hues indicated by BCS, and choose the 
most favorable. 

An analysis of this sort indicates that th£ true ground state contains 



We have explicitly displayed color (a,/3), flavor (i, j), and spinor (a, 6) indices. 
The A-index is summed and therefore links color and flavor. We have used 
a two-component spinor notation; note that properly the right-helicity fields 
should involve dotted spinors. The important information conveyed by the 
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spinors is that the condensation does not violate rotational invariance. The 
relative minus sign between left-helicity and right-helicity condensates signifies 
that the ground state is a scalar, rather than a pseudoscalar, so that parity 
is unbroken. The magnitude of the condensate depends on the magnitude of 
the 3-momentum p. In weak coupling, it is largest for |p| near the nominal 
(free-particle) Fermi surface. 

Many different treatments have shown thaJjpa_c£ffl|depsate of the form (|l|) is 
the dominant condensate in three-flavor QCDQ'Bll3llil^t3 The essential physical 
argument that favors this pattern is that, by leaving the maximal unbroken 
symmetry group, this pattern allows quarks of all three colors and all three 
flavors to pairlJ Less symmetric condensates in which not all ouarks pair do 
not lower the free energy as much as the "maximal" choice (|l|)ll3't3 

In reality, condensation in the color 3 channel (|^) induces a small but 
nonzero condensate|-in the color 6 channel even if the interaction is repul- 
sive in this_channel0 because this additional condensation breaks no further 
symmetriesllj This means that the right hand side of (|l|) is slightly more com- 
plicated and should, in fact, be written in terms of two gap parameters ki and 
K2, as Ki(p^)5"(5^ -I- K2(p^)S^S^- The pure color 3 condensate displayed in 
has K2 = —Ki. Using (|l|) is a good approximation because the induced color 6 
condensate is much smaller tharLthfi-dominant color 3 condensate mandated 
by the attraction in this channelQEJO 

We can now explain the term "color-flavor locking" . Writing e^^^eabA — 
^a^b ~^b^a^ see that the condensates (0) involve Kronecker delta functions 
that link color and flavor indices. These condensates transform nontrivially 
under separate color and flavor transformations. Neither color transformations 
nor flavor transformations, separately, are valid symmetries of the ground state. 
However, the delta functions do remain invariant if we simultaneously rotate 
both color and flavor. Thus these symmetries are locked together. 

Pairing which locks two previously unrelated symmetries is not a new 
phenomenon. The condensate for the color-flavor locked state of QCD perhaps 
most directly resembles that found in the B phase of superfluid helium 3. 
In an ordinary nonrelativistic system, orbital and spin rotations are separate 
symmetries. In the B phase, however, a condensate of pairs of atoms formed 
at the Fermi surface transforms nontrivially under spin and orbital rotations 
but is invariant under simultaneous rotations of both: it locks spin and orbital 
rotations. Another analogy is to the physics of chiral symmetry breaking in 
the QCD vacuum. There, the formation of a condensate of left-handed quarks 
and right-handed antiquarks locks the S'C/(3)l and SU{?i)ji flavor symmetries. 

All the main qualitative properties of the color-flavor locked state are direct 
consequences of (^, properly interpreted. Elucidation of these consequences 
will occupy us for the remainder of the Section. 



The Condensed Matter Physics of QCD 



9 



2.2 Symmetry Breaking 

An aspect of (|^) that might appear troubling at first sight is its lack of gauge 
invariance. This actually turns out to be a profound advantage. 

TherCpare powerful general arguments that local gauge invariance cannot 
be brokenEJ Indeed, local gauge invariance is really a tautology, stating the 
redundancy of variables. Yet its "breaking" is central to two of the most 
successful theories in physics, namely BCS superconductivity theory and the 
standard model of electroweak interactions. In BCS theory we postulate a 
nonzero vacuum expectation value for the Cooper pair field, which is elec- 
trically charged. In the electroweak standard model we postulate a nonzero 
value for the Higgs field, which violates both the weak isospin SU{2) and the 
hypercharge U{1) gauge symmetries. 

In each case, we should interpret the condensate as follows. We are working 
in a gauge theory at weak coupling. It is then very convenient to fix a gauge, 
because after we have done so — but not before! — the gauge potentials in 
which we perturb will make only small fluctuations around zero. Of course at 
the end of any calculation we must restore the gauge symmetry, by averaging 
over the gauge fixing parameters. Only gauge-invariant results will survive 
this averaging. However, in the intermediate steps, within a fixed gauge, one 
can capture important correlations that characterize the ground state by spec- 
ifying the existence of nonzero condensates relative to the ambient gauge. In 
superconductivity, the essence of the physics is the correlation in the fermionic 
wave function which describes the Cooper pairs, and the resulting modification 
of the dispersion relations which describe the excitation spectrum. In partic- 
ular, the gap in the spectrum of fermionic excitations at the Fermi surface is 
a gauge invariant quantity. Describing this physics within a fixed gauge as a 
condensate which "breaks" the gauge symmetry is a convenient fiction. 

In the standard electroweak model one employs a nonzero vacuum expec- 
tation value for a Higgs doublet field {(j>°') = vS^, which is not gauge invariant. 
One might be tempted to use the magnitude of its absolute square, which is 
gauge invariant, as an order parameter to characterize the symmetry breaking. 
However, (0^0) never vanishes, whether or not any symmetry is broken. In 
fact there can be no gauge invariant order parameter for the electroweak phase 
transition, since it has long been known that one can, by allowing the SU(2) 
gauge couplings to become large, go over intp-a "confined" regime, encounter- 
ing no sharp phase transition along the wayuJ The absence of massless gauge 
bosons and of long-range forces is the essence of the Meissner-Anderson-Higgs 
effect — and it is also the essence of confinement! 

So evidently, when used with care, the notion of spontaneous gauge sym- 
metry breaking can be an extremely convenient fiction. In particular, by forg- 
ing a connection with superconductivity and condensate formation, it brings 
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the universality class of confinement down to earth, and makes it accessible 
to weak coupling methods. These condensates need not break any true (i.e. 
global) symmetries. If a global symmetry is broken, some combination of the 
condensates themselves is a gauge invariant physical observable, and not just 
a convenient fiction. 

With this discussion in mind, let us consider the consequences of (|l|) for 
symmetriesLl (See also Ref. 16 in which this ordering was considered at zero 
density.) 

The exact microscopic symmetries of QCD with three massless flavors are 

^^microscopic — 5C/(3)color X SU{i)L X SU{i)R X U{1)b , (2) 

where the first factor is the (local) vector color gauge symmetry, the second and 
third factors are global chiral flavor symmetries, and the fourth factor is baryon 
number. (For the present, we set the electromagnetic coupling constant to zero, 
and therefore neglect the fact that the U{1) subgroup of SU{3)l x SU{3)ii 
which describes electromagnetism is in fact a local symmetry.) 

The color-flavor locked phase (||) features two condensates, one involving 
left-handed quarks alone and one involving right-handed quarks alone. The 
former locks SU{3)l flavor rotations to SU {3)coIot- the condensate is not sym- 
metric under either alone, but is symmetric under simultaneous iS'C/(3)L+coior 
rotations. Similarly, the condensate involving right-handed quarks alone locks 
SU{3)ii flavor rotations to SU{3)coior- As a consequence, of all the symmetries 
in Gmicroscopic, Only the subgroup 

GcFL = <S'C/(3)color+L+i?. X Z2 (3) 

leaves the correlated ground state (|^) invariant. The color and chiral flavor 
symmetries are broken, by color-flavor locking, down to the (global) vector 
"diagonal" symmetry, that makes equal transformations in all three sectors - 
color, left-handed flavor, and right-handed flavor. Baryon number symmetry 
is broken down to a discrete Z2 symmetry under which all quark fields are 
multiplied by —1. Even though the condensates in (|l]) do not appear to lock 
SU{3)l to SU (3)_R, they manage to do so by locking both to SU{3)co\or- Color- 
flavor locking, therefore, provides a mechanism by which chiral symmetry can 
be broken. 

The spontaneous breaking of chiral symmetry is a familiar phenomenon 
in zero-density QCD. Here, at high density, it occurs by a rather different 
mechanism. In zero-density QCD chiral symmetry breaking is due to con- 
densation of left-handed quarks with right-handed antiquarks. The pairing of 
opposite helicities, of course, breaks chiral symmetry. Here we have only pair- 
ing of left-handed quarks with left-handed quarks, and right-handed quarks 
with right-handed quarks. Nevertheless chiral symmetry is broken indirectly, 
as we have described. 
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Now we must mention explicitly a general formal consequence of symmetry 
breaking, that we glided over earlier. That is, the form of the condensate in 
is not unique. There are states with equally good energy where the correlated 
ground state is subjected to the action of any transformation in Gmicroscopic- 
This action will in general (except for elements of Gcfl) alter the form of 
the condensate. So there is a manifold of distinct vacua associated with the 
broken symmetries. By interpolating among these vacua, with small space- 
time gradients, one produces low-energy excitations around a fixed, reference 
ground state such as the one specified by (|l|). 

The significance of such states is familiar from other contexts. When 
the broken symmetry is a global symmetry, they lead to Nambu-Goldstone 
particles. When the broken symmetry is a local gauge symmetry, the Meissner- 
Anderson-Higgs mechanism is triggered, and the would-be Nambu-Goldstone 
particles become the longitudinal parts of massive gauge fields. 

Here the local color symmetry has been broken completely. Thus all the 
gluons acquire mass. This result has a major positive consequence for the 
logical status of our analysis. It removes the possibility of infrared divergences 
associated with gluon exchange. Similarly, the existence of an energy gap for 
all the quarks removes the other potential source of infrared divergences, from 
integration over low-energy excitations around the Fermi surface. Altogether, 
then, weak coupling perturbation theory around the correct, condensed ground 
state is free of the difficulties that appeared around the naive ground state. 

The spontaneous violation of baryon number is perhaps less familiar to a 
particle physicist. On first hearing, one might think this is a dramatic or even 
catastrophic prediction, since we know that baryonic matter in the Universe 
is stable over very long periods of time. That is obviously too naive an inter- 
pretation, however, since in the theory of ordinary superconductors we deal 
with electron pairing without worrying over the violation of lepton number, 
and in the theory of helium superfluids we deal with condensates of atoms or 
diatoms, which formally violate both baryon and lepton number. Finally, and 
(as we shall see below) most directly to the point, ordinary nuclear matter is a 
superfluid in which nucleon-nucleon pairing violates baryon number symmetry. 

In all these cases, if we are dealing with a finite sample of the superconduc- 
tor or superfluid, there is no true violation of the conservation laws. Indeed, 
we may draw a surface surrounding the sample, and apply Gauss' law to the 
equation of current conservation (modified, within the sample, to include the 
condensate) in the usual way to see that changes in the bulk quantum num- 
bers are accompanied by appropriate compensating fluxes. The correct inter- 
pretation of the formal "violation" of these symmetries is that there can be 
large fluctuations and easy transport of the corresponding quantum numbers 
within the sample. These are precisely the phenomena of superconductivity 
and superfluidity. The mathematical connection between broken symmetry 
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and super-transport is quite direct: As we have discussed, the symmetry break- 
ing order parameter can point in any one of a manifold of "directions" ; the 
supercurrents are carried by spatial variation in the "direction" of the conden- 
sate; that is, the phenomenon of superfluidity is a direct manifestation of the 
Nambu-Goldstone mode associated with the spontaneous breaking of a global 
symmetry; similarly, superconductivity is a direct manifestation of the mode 
which becomes the longitudinal component of a massive gauge boson and is 
thus a direct manifestation of the breaking of a gauge symmetry. 

As already mentioned, the standard Higgs mechanism in the electroweak 
sector of the standard model has no gauge-invariant order parameter. With 
color-flavor locking the situation is more fortunate, because global as well as 
gauge symmetries are broken. Physically, this implies that there are sharp 
differences between the color-flavor locked phase and the quark-gluon plasma 
phase (in which all symmetries of the QCD lagrangian are unbroken), so that 
any passage between them must be marked by one or more phase transitions. 
In fact, it is a simple matter to abstract gauge invariant order parameters, 
which have a strict meaning valid at any coupling, from our primary, gauge 
variant condensate at weak coupling. For instance, to form a gauge invariant 
order parameter capturing chiral symmetry breaking we may take the product 
of the left-handed version of (|l|) with the right-handed version and saturate 
the color indices, to obtain 

{rui^lji^RjRp) - {ru^l,){i'RjRp) - A26,,„,e'='™ . (4) 

Likewise we can take a product of three copies of the condensate and sat- 
urate the color indices, to obtain a gauge invariant order parameter for the 
baryon-number violating superfluid order parameter. These secondary order 
parameters will survive gauge unfixing unscathed. Unlike the primary conden- 
sate, from which they were derived, they are more than convenient fictions. 

If we turn on a common mass for all the quarks, the chiral SU (3) l x SU (3)_r 
flavor symmetry of Gmicroscopic will be reduced to the diagonal SU(i)L+R- If 
we turn on unequal masses, the symmetry will be even less. In any case, 
however, the ^7(1) of baryon number a good microscopic symmetry, and the 
corresponding six-quark order parameter remains a strict signature of the color- 
flavor locked phase, distinguishing it both from the quark-gluon plasma phase, 
and from some other states of quark matter we shall encounter in Section 3. 

As it stands the order parameter (Q) is not quite the usual one, but roughly 
speaking its square. It leaves invariant an additional Z2, under which only 
the left-handed quark fields change sign. Actually this Z2 is not a legitimate 
symmetry of the full theory, but suffers from an anomaly. So we might expect 
that the usual chiral order parameter is induced by the anomalous interactions 
that violate the axial baryon number symmetry of the classical Lagrangian. To 
put this another way, because axial baryon number is not a symmetry of QCD, 
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once chiral symmetry is broken by color-flavor locking there is no symmetry 
argument precluding the existence of an ordinary chiral condensate. Indeed, 
instanton elTects do induce a nonzero {'iPriPl) because the instanton-induced 
interaction is a six-fermion operator which can be written as a product of '^RipL 
and the operator in which already has a nonzero expectation valueP but 
this turns out to be a small effectaM 

At weak coupling, we can be more specific about these matters. The 
most important interactions near the Fermi surface, quantitatively, arise from 
gluon exchange. These are responsible for the primary condensation. The in- 
stanton interaction is much less important asymptotically because the gauge 
fields which make up the instantons themselves are screened, the effects of 
instantons are intrinsically smaller for more energetic quarks, and because 
the instanton-induced interaction involves six fermion fields, and hence (one 
can show) becomes irrelevant upon renormalization toward the Fermi surface. 
The instanton interaction is qualitatively important, however, because it rep- 
resents the leading contribution to axial baryon number violation. It is only 
such U{1)a violating interactions that remove the degeneracy among states 
with different relative phases between the left- and right-handed condensates 
in (|l|). In the absence of intrinsic U{1)a breaking, the spontaneous violation of 
this symmetry in the color-flavor locked phase would be accompanied by the 
existence of a pseudoscalar SU{3)coIot+l+r singlet Nambu-Goldstone bosons. 
Since the intrinsic violation of this symmetry is parametrically small, the cor- 
responding boson will not be strictly massless, but only very light. Quantum 
fluctuations in this light 77'-field, among other things, will keep the conven- 
tional chiral symmetry breaking order parameter small compared to at 
high density. 



2.3 Elementary Excitations 



The physics J^l£p^^j&,eJpita.tkjpS||^ much 
recent workByy y yy'yHH'y'Hkilfea'BEa There 
are three sorts of elementary excitations. They are the modes produced di- 
rectly by the fundamental quark and gluon fields, and the collective modes 
associated with spontaneous symmetry breaking. These modes can be classi- 
fied under the unbroken SU (3) x Z2 symmetry, and the unbroken rotation and 
parity symmetries. 

The quark fields of course produce spin 1/2 fcrmions. Some of these are 
true long-lived quasiparticles, since there are no lighter states of half-integer 
spin that they might decay into. With the conventions we have been using, as 
expressed in ([l]), the quark fields are triplets and antitriplets under color and 
flavor, respectively. Thus they decompose into an octet and a singlet under 
the diagonal S'?/(3)coior+L+_R- There is an energy gap for production of pairs 



14 



Handbook of QCD / Volume 3 



above the ground state. More precisely, tbere are two gaps: a smaller one for 
the octet, and a larger one for the singlet Q The dispersion relations describing 
these ferrnionic_aiiaaiparticle excitations in the CFL phase have been described 
in some detailEJcjO 

The gluon fields produce an SU{2>)co\or+L+R octet of spin 1 bosons. As 
previously mentioned, they acquire a common mass by the Meissner-Anderson- 
Higgs mechanism. The quantitative expressions for the maaaea^pf these vec- 
tor mesons which havejpeen computed at weak coupling o'ESEaa and in an 
instanton-liquid modelED 

The fermionic excitations have a gap; the vector mesons have mass; but, 
the Nambu-Goldstone bosons are massless. These bosons form a pseudoscalar 
octet associated with chiral symmetry breaking, and a scalar singlet associated 
with baryon number superfluidity. The octet, but not the singlet, will be lifted 
from zero mass if the quarks are massive. Finally there is the parametrically 
light, but never strictly massless, pseudoscalar singlet associated with U{1)a 
breaking. 

The Nambu-Goldstone bosons arising from chiral symmetry breaking in 
the CFL phase are Fermi surface excitations in which the orientation of the 
left-handed and right-handed diquark condensates oscillate out of phase in 
flavor spacCj— .TJae-cffective field theory describing these oscillations has been 
constructed HlUcj Up to two derivatives, it is given by 

r.ff = :^Tr (aoSaoS^ + vld^m'Y.'^) - c (detMTr(M-iE) + h.c.) . (5) 

The Nambu-Goldstone boson field matrix S is a color singlet and transforms 
under SU{3)l x SU{3)r as E — > Ul^uI^ as usual. M — diag(TO„, md, TOs) is 
the quark mass matrix. The construction of E from rotations of the CFL con- 
densates can be found in Refs. 19,21: one first finds the 19 putative Nambu- 
Goldstone bosons expected when Gmicroscopic is broken to Gcfl; one then 
identifies 8 of these which become the longitudinal parts of massive vector 
bosons; the remaining ten are the octet described by (^, the singlet 77' and 
the superfiuid mode. See Refs. 19,21 for the singlet terms in the effective La- 
grangian. The_higher derivative terms in the effective Lagrangian have also 
been analyzedE^ It has also been suggested that once higher derivative inter- 
actions are introduced, the effective theory may support Skyrmions, in which 
the Nambu-Goldstone boson field configuration has rLOUzero winding numberlia 
These solitons have energies comparable to the gapp3 and are an alternative 
description of the gapped fermionic excitations in the CFL phase, in the same 
sense that baryons can alternatively be described as Skyrmions of the vacuum 
pion field. 

The masses of the pseudoscalar mesons which are the pseudo-Nambu- 
Goldstone bosons associated with chiral symmetry breaking can be obtained 
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from £off of (|5|)Lil For example 
2c 

J IT TV 



m^± = -^ms{mu + rud) , m^± = -7jmd(m„ + m^) . (6) 



Thus, the kaon is Ughter than the pion, by a factor of rrid/imu + md)^ Note 
that the effective Lagrangian is quadratic in AI. This arises because Leff ^' 
spects the Z2 symmetry under which only the left-handed quarks change signl 
As we discussed in the previous section, this is almost a symmetry of thfi CFL 
phase: it would be a symmetry if instanton effects could be neglectedEl How- 
ever, instanton effects genemte a nonzero, but small, ordinary (ipRipL) conden- 
sate, which breaks the Z2 0'lI'c3 and results in a contribution to-ihe meson 
which is linear in M and which may be numericalhj-aignificantcfLThe induced 
{'4'R'^l) was significantly overestimated in Refs. lI'E£|, howeverE3 Making the 
already small {iPr'4'l) even smaller will significantly reduce the Z2-violating 
contributions to the meson masses. 

If we were describing pions in vacuum, or pions in nuclear matter, the 
only way to obtain the cocfhcients in the effective theory would be to measure 
them in an experiment or, if possible, to calculate them on the lattice. Indeed 
in any theory with strong interactions, the purpose of writing an effective 
theory for the low energy degrees of freedom is to express the predictions 
for many low energy processes in terms of a few parameters, which must be 
obtained from experiment. In the color-flavor locked phase, however, the full 
theory is weakly coupled at asymptotically high densities. In this regime, 
therefore, the coefficients and c are calculable from first principles using 

ijpaV rniinling TTipfVinrlal Up to possible logarithmic corrections, the result 



21-81og2 1 3A^ 

Sdtt^ 3 27r^ 



The electromagnetiSH and nonzero temperaturS corrections to these quan- 
tities have also been calculated. 

Quantitatively, (see Section 4 for a discussion of estimates of A) we esti- 
mate that the lightest pseudoscalar meson, the kaon, has a mass in the range 
of 5 to 20 MeV at /i = 400 MeV, and becomes lighter still at higher densities. 
There are two reasons why the Nambu-Goldstone bosons are so much lighter in 
the CFL phase than in the vacuum. First, their mass^ is proportional to m'^^^^.^^ 
rather than to JTiquark, as at zero density. In addition, there is a further sup- 
pression by a factor of A/n, which arises because the Nambu-Goldstone bosons 
are collective excitations of the condensates formed from particle-particle and 
hole-hole pairs near the Fermi surface, whereas the quark mass term connects 
particles with antiparticles, far from the Fermi surfaceE 
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In QCD with unequal quark masses, at very high densities the CFL phase 
is much as we have described it, except that the gaps associated with (its), 
(ds) and {ud) pairing wiU difFer-shghtly, the Fermi momenta for the different 
quark flavors may differ slightlyP3 and there may therefore be a small electron 
chemical potential /ig. (We discuss lower densities, where the differences be- 
tween Fermi momenta become comparable to A, in Section 3.) Because the 
strange quarks are less numerous due to their greater mass, it would seem that 
the requirement of overall charge neutrality can only be satisfied if /ie > nie 
and there is a nonzero density of electrons. Schafer has recently noted, how- 
ever, that if /ie > mx-, charge neutrality will be achieved by the formation 
of a condensate oi-K^ bosons, rather than by the formation of a small Fermi 
sphere of electronsE3 The CFL condensate rotates in the direction, relative 
to the direction favored by the quark mass term in the effective Lagrangian. 
This costs energy, but because it introduces a negative electric charge it is 
free-energetically favored, if /ie > to_r-- . 

The formation of a kaon condensate in the CFL phase changes none of its 
symmetries. In this sense, it is a less dramatic effect than the formation of a 
kaon condensate in nuclear matter made of neutrons and protons only:E3 there, 
kaon condensation breaks U{l)s- This symmetry is already broken in the CFL 
phase. Kaon condensation in the CFL phase is more akin to kaon condensation 
in hypernuclear matter made up of equal measures of all the octet baryons, in 
which U{l)s is already broken by hyperon-hyperon pairing. We shall elaborate 
much further on this connection in Section 2.5. 

2.4 The Modification of Electromagnetism 

It is physically significant, and proves extremely instructive, to consider the 
effect of color-flavor locking on the electromagnetic properties of high-density 
hadronic matter. 

To do this, we consider coupling in the appropriate additional C/(1)em 
gauge field A^, representing the photon. This couples to u,d,s quarks with 
strength |e, — ^e, — ^e, respectively. Evidently this J7(1)em symmetry is bro- 
ken by the condensate (^, through the terms pairing differently-charged quarks. 
Were this the complete story, the color-flavor locked phase would be an elec- 
tromagnetic superconductor. The truth is far different, however. 

The situation is analogous to what occurs in the electroweak sector of the 
standard model. There, the Higgs field condensate breaks both the original 
weak SU{2) and the hypercharge U{1). However, one linear combination of 
these symmetries leaves the condensate invariant, and remains a valid gauge 
symmetry of the ground state. Indeed, this is how we identify electromagnetism 
within the standard model. 

Here we must similarly consider the possibility of mixing broken color 
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SU{?>) and (original) electromagnetic U{1)-e,m generators to find a valid residual 
symmetry. Indeed, we should expect this to occur, by the following argument. 
In QCD with three flavors, J7(1)em is a subgroup of SU{i)Lj^f!. When we 
neglected electromagnetism, we found that in the color-flavor locked phase 
S'J7(3)i4.fl is broken but S'f/(3)coior+L+R is an unbroken global symmetry. We 
therefore expect that gauging a U{1) subgroup of SU{2>)l+r must correspond, 
in the CFL phase, to gauging a U{1) subgroup of the unbroken S't/(3)coior+L+R- 
Once we are alerted to this possibility, it is not difficult to identify Ihe 
appropriate combination of the photon and gluons which remains unbrokenDEJ 
The unbroken U{1) is generated by 

Q^Q + TjTs (8) 

with 77 — l/-\/3- Q is the conventional electromagnetic charge generator and 
Tg is associated with one of the gluons. In the representation of the quarks, 

Q — diag(|, — i, — i) in flavor u, d, s space, 

(9) 

Tg — diag(— 2, 1, 1) in color r, 5, 6 space. 

As is conventional, we have taken tr(T8r8) = 2. By construction, the Q-charges 
of all the Cooper pairs in the condensate (|l]) are zero. (For example, with these 
conventions, red up quarks pair only with green down or blue strange quarks, 
and both these pairs have Q = in sum.) The condensate is Q-neutral, the 
U{1) symmetry generated by Q is unbroken, and the associated Q-photon will 
remain massless. To see exactly which gauge field remains unbroken, look at 
the covariant derivative of the condensate, 

D^MaQ^^) - (a^ + eA,,Q + gGlTs) (q^q^) . (10) 

The kinetic term |-D(?a?f)P will give a mass to one gauge field 

„ —rieA,, + oG?, 

= ' ^\ / - - sin^o^p + cosaoG^ , (11) 

while the orthogonal linear combination 



AQ = ^ ' ; = cosaoA^ +sinaoG« (12) 

is the Q-photon which remains massless. That is, B'^ satisfies the ordinary 
Maxwell equations while experiences a Meissner effect. The denominators 
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arise from keeping the gauge field kinetic terms correctly normalized, and we 
have defined the angle ao which specifies the unbroken U{1) via 

cosao — ^ = ■ (13) 

The mixing angle ag is the analogue of the Weinberg angle in electroweak 
theory, in which the presence of the Higgs condensate causes the and the 
third SU{2)y/ gauge boson to mix to form the photon, A^, and the massive Z 
boson. At accessible densities the gluons are strongly coupled (g^/(47r) 1), 
and of course the photons are weakly coupled (e^/(47r) « 1/137), so Uq ~ r]e/g 
is small. The "rotated photon" consists mostly of the usual photon, with only 
a small admixture of the Tg gluon. 

Let us now consider the charges of all the elementary excitations which 
we enumerated in Section 2.3. For reference, the electron couples to with 
charge 

(14) 



which is less than e because the electron couples only to the component of 
. Now in computing the Q-charge of the quark with color and flavor indices 
a, a we must take the appropriate combination from 

e(-f ff' iff' jg) + g(ie^ -\e^ -\) 

One readily perceives that the possible values are 0, ±e. Thus, in units of the 
electron charge, the quarks carry integer Q-charge! Quite remarkably, high- 
density QCD realizes a mathematically consistent gauge theory version of the 
old vision of Han and Nambu: the physical quark excitations have integer 
electric charges that depend on an internal color quantum number! 

Similarly, the gluons all have Q-charges 0,±e. Indeed, they have the Q- 
charges one would expect for an octet of massive vector bosons. The Nambu- 
Goldstone bosons arising from the breaking of chiral symmetry, of course, have 
the same charge assignments as the familiar tt, K and 77 octet of pseudoscalars. 
The baryon superfluid mode is Q-neutral. In the color-flavor locked phase, we 
conclude, all the elementary excitations are integrally charged.]^ This is a 
classic aspect of confinement, here embodied in a controlled, weak-coupling 
framework. 



"We shall see in Section 3 that in two-flavor QCD, in which color-flavor locking does not 
occur, the color superconducting condensate which forms also leaves a Q-photon massless. 
The only difference relative to the CFL phase is that r) = — l/2v^. (However, the Q-charges 
of the excitations are not all integral in this theory.) 
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It is fun to consider how a chunk of our color-flavor locked material would 
look. If the quarks were truly massless, then so would be some charged Nambu- 
Goldstone bosons, and one might expect a rather unusual "bosonic metal" , in 
which the low-energy electromagnetic response is dominated by these modes. 
As we mentioned in Section 2.3, turning on quark masses gives small masses to 
the Nambu-Goldstone bosons. Electromagnetic radiative corrections further 
lift the masses of the charged Nambu-Goldstone bosons. As a result, the 
color-flavor locked material becomes a transparent insulator, with no charged 
excitations at zero temperature. Altogether, a chunk of color-flavor locked 
material would resemble a diamond: an ordinary light wave incident upon it 
would be partially reflected, but some fraction would be admitted as a Q- 
light wave; the Q-light would travel through the transparent "diamond" , and 
would partially emerge and be partially internally reflected on the far side. 
Turning on small but unequal quark masses mars the diamond somewhat: Ea 
maintaining overall electric charge neutrality in this casa-sYill require either a 
nonzero electron density or a condensate of charged kaonspO and in either case, 
Q-light will have charged excitations off which to scatter. 

Although a quantitative calculation of light reflecting off the facets of a 
CFL-diamond has not yet been done, the effect of a chunk of color supercon- 
ducting quark matter (whether in the CFL phase or in the less symmetric 
phase in which only up and down nuarks pair) on a static magnetic fleld has 
been described in complete details Some fraction of an externally applied 
ordinary magnetic field penetrates the superconductor in the form of a Q- 
magnetic field, while some fraction of the ordinary magnetic field is expelled 
by the Meissner effect. The fraction of the field which is expelled depends both 
on ao and on the shape of the chunk color superconducting quark matter, but 
it is small when ag is small, as in nature. Most of the flux is admitted, as Q- 
flux. This Q-magnetic field satisfies Maxwell's equations and is not restricted 
to flux tubes. 



2.5 Quark-Hadron Continuity 

The universal features of the color-flavor locked state: confinement, chiral sym- 
metry breaking leaving a vector SU (3) unbroken, and baryon number super- 
fluidity, are exactly what one expects to find in nuclear matter in three-flavor 
QCDeI Perhaps this is not immediately obvious in the case of baryon number 
superfluidity, but let us recall that pairing phenomena, which would go over 
into neutron superfluidity and proton superconductivity in nuclear matter, are 
very well established in ordinary nuclei. In three-flavor QCD, there are good 
reasoned to think that the pairing interaction in the flavor singlet dibaryon 
channel (the so-called 7J-dibaryon channel) would be quite attractive in three- 
flavor QCD, and support a robust baryon number superfluidity. Thus, the 
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symmetries of the color-flavor locked phase are precisely those of nuclear mat- 
ter in three- flavor QCD, perhaps better referred to as hypernuclear matterH 

Furthermore, there is an uncanny resemblance between the low-lying spec- 
trum computed from first principles for QCD at asymptotically high density, 
and what one expects to find in hypernuclear matter, in a world with three 
degenerate quark flavors. It is hard to resist the inference that in this theorv. 
there need be no phase transistion between nuclear density and high densityB 
There need be no sharp boundary between hypernuclear matter, where micro- 
scopic caculations are difficult but the convenient degrees of freedom are "obvi- 
ously" hadrons, and the asymptotic high-density phase, where weak-coupling 
(but non-perturbative) calculations are possible, and the right degrees of free- 
dom are elementary quarks and gluons, together with the. collective Nambu- 
Goldstone modes. We call this quark-hadron continuityH Perhaps the least 
surprising aspect of this, by now, is the continuity between the pseudoscalar 
mesons at nuclear density and those at asymptotically high densities, since in 
both regimes these are present as a result of the breaking of the same symme- 
try. It might seem more shocking that a quark can go over continuously into, 
or "be" , a baryon, since baryons are supposed to contain three quarks, but 
remember that in the color-flavor locked phase the quarks are immersed in a 
diquark condensate, and so a one-quark excitation can pick two quarks up from 
(or lose to quarks to) the condensate at will. The difference between one and 
three is negotiable. What about the gluons? Within the color-flavor locked 
phase, similarly, they are quite literally the physical vector mesons. They are 
massive, as we have discussed, and have the right quantum numbers. Thus the 
original vision of Yang and Mills - who proposed non-abelian gauge theory as 
a model of p mesons - is here embodied. 

Note that the hypothesis of continuity between hypernuclear and dense 
quark matter certainly does not preclude quantitative change. Far from it. 
The dispersion relation for a fermion — whether a quark in the CFL phase or 
a baryon in the hypernuclear phase — is characterized by a gap at the Fermi 
surface and by a gap at zero momentum, i.e. a mass. As a function of increasing 
density, gaps at the hyperon Fermi surfaces due to hyperon-hyperon pairing 
evolve continuously to become the gaps at, the quark Fermi surfaces which 
characterize the color-flavor locked phaseEJ During this evolution, the gaps 
are thought to increase significantly. In contrast, the gaps at zero momentum 
decrease dramatically with increasing density as they evolve from being of 
order the hyperon masses in hypernuclear matter to being of order the current 
quark masses at asymptotically high densities. 

Note that in order for quark-hadron contirmity to be realized, U{\)em 
must not be broken by hyperon-hyperon pairingEj At every point during the 
evolution of the theory as a function of increasing density, there will be an 
unbroken U(V) and a massless Q-photon and the excitations will be integer 
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charged. As the density is increased, however, the definition of the Q-photon 
in terms of the vacuum photon and gluon fields rotates. The ordinary photon 
rotates to become the Q-photon of the CFL phase, defined in (jl^). Turning 
to the massive vector bosons, in hypernuclear matter there is both an octet 
and a singlet. The singlet must become much heavier than the octet as a 
function of increasing density, since in the low energy description of the color- 
fiavor locked phase one finds the octet alone. We see that if quark-hadron 
continuity is realized in QCD with three degenerate quarks, it requires various 
quantititative (but continuous) changes. What is remarkable is that it is even 
possible to imagine watching all the physical excitations of the theory evolving 
continuously as one dials the density up and goes from a strongly coupled 
hadronic world to a weak-coupling world of quarks and gluons. 

If the quarks are massless, the Nambu-Goldstone bosons are massless in 
both hypernuclear and CFL quark matter, and in between. Once nonde- 
generate quark masses are introduced, however, the evolution of the Nambu- 
Goldstone masses as a function of increasing density becomes more intricate, 
as the kaon must go from being heavier than the pion to being lighter. 

Note, finally, that the whole story becomes^further complicated once the 
strange quark is made as heavy as in natureli^O Although the color-flavor 
locked phase is certainly obtained at asymptotically densities, where quark 
masses are neglectable, the nuclear matter phase, made of neutrons and pro- 
tons only, is not continuously connectable with the color-flavor locked phase. If 
quark-hadron continuity is to be realized in the phase diagram of nature, what 
must happen is that, as a function of increasing density, one first goes from nu- 
clear matter to hypernuclear matter, with sufficiently high density that all the 
hyperons have similar Fermi surfaces. This first stage must involve phase tran- 
sitions, as the symmetries of hypernuclear matter differ from those of ordinary 
nuclear matter. Then, as the density is increased further, the hypernuclear 
matter may evolve continuosly to become CFL quark matter, with pairing 
among hyperons becoming CFL pairing among quarks. 

We now have a description of the properties of the CFL phase and its 
excitations, in which much can be described quantitatively if the value of the 
gap A is known. We describe estimates of A in Section 4. First, however, we 
give a full description of the less symmetric variants of color superconductivity 
which arise in QCD with Nf ^ 3. Already, however, in our idealized world 
(in which we either have three degenerate quarks or such high densities that 
the quark mass differences can be neglected) let us pause to marvel at our 
theoretical good fortune. The color-flavor locked phase is a concrete realization 
of the idea of complementarity: the same phase of a gauge theory can be 
described simultaneously as one in which the^gauge symmetry is spontaneously 
broken and as one in which color is confinedllj This means that it provides us 
with a weak-coupling laboratory within which we can study a confined phase 
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from first principles at weak coupling. It is furthermore a phase of QCD 
wherein the physics of chiral symmetry breaking — indeed all the parameters 
of the chiral effective Lagrangian and all known or conjectured phenomena of 
the pseudoscalar meson sector, including kaon condensation — are amenable 
to controlled, weak-coupling calculation. 



3 2SC and Other Variants 

3.1 Two Flavors 

In the previous Section, we have described quark matter in QCD with three 
degenerate flavors of light quarks. Nature is less symmetric, and in order to 
bracket nature we now describe the color superconducting phase in QCD with 
two flavors of light quarks. Pairs of quarks cannot be color singlets, and in 



QCD with two flavors o: 
tractive) color 3 channel 



massless quarks the Cooper pairs form in the (at- 
'Bi'l The resulting (eapa^^^^f^'j) condensate picks 
a color direction (in this case the 3 or blue direction), creates a gap A at the 
Fermi surfaces of quarks with the other two out of three colors (red and green) , 
and breaks SU{3)co\ot to an S'[/(2)coior subgroup, giving mass to five of the 
gluons by the Anderson-Higgs mechanism. The masses of these vector bosons 
have beeiLjComputed in the weak-coupling theory valid at asymptotically high 
densitiesfj and in an instanton liquid modelEll Axial color is not a symme- 
try of the QCD action, but at asymptotically high densities where the QCD 
coupling g is weak, explicit axial color breaking is also weak. As a result, the 
pseudoscalar excitations of the condensate which would be Goldstone bosons 
arising from axial-S'f7(3)coior to axial-S'f7(2)coior breaking if g were zero may 
be rather lightca 

In QCD with two flavors, the Cooper pairs are ud — du flavor singlets 
and the global flavor symmetry SU{2)l x SU{2)ii is intact. There is also an 
unbroken global symmetry which plays the role of U{1)b- Thus, no global 
symmetries are broken in this 2SC phase. This means that as a function of 
increasing density, there must be a phase transition between the hadronic and 
2SC phaspji; at which chiral symmetry is restored. This phase transition is first 
ordetlo'EfrLJ since it |ixi¥Qlves a competition between chiral condensation and 
diquark condensationc3'tj 

Because no global symmetries are broken, there are no light scalar degrees 
of freedom in the 2SC phase. The 2SC phase is not a superfluid. The poten- 
tially light degrees of freedom are the gapless blue quarks, and the gauge bosons 
associated with unbroken gauge symmetries. The gluons associated with the 
unbroken S'[/(2)coior will exhibit strong dynamics on long enough length scales; 
this aspect of the infrared physics of the 2SC phase is not under perturbative 
control. Thus, whereas in the CFL phase we expect that any physical quan- 
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tity can be obtained from a controlled weak-coupling calculation at sufficiently 
high density, this claim cannot be made for the 2SC phase. Note that the 
gapless blue quarks are neutral under S'J7(2)coior- As in the CFL phase, there 
is a massless abelian gauge boson, formed from a suitable linear combination 
of the photon and one of the eight original gluons. The two blue quarks have 
charges and 1 under this unbroken, but rotated, Q-electromagnetism. 

We expect that the blue quarks, left ungapped by the primary 2SC con- 
densate, will find some way to pair in a higher angular momentum channel. 
Indeed, the instanton interaction pairs these quarks in a J = 1 condensate 
which breaks rotational invariance. Early work suggested that the associated 
gap is of order keVpl but this estimate should be revisited. 

It is interesting that both the 2SC and CFL phases satisfy anomaly match- 
ing constraints, even though it is not yet completely clear whetheiLthis must be 
the case when Lorentz invariance is broken by a nonzero depsitycS It is not yet 
clear how high density QCD with larger numbers of flavorsP^ which we discuss 
below, satisfies anomaly matching constraints. Also, anomaly matching in the 
2SC phase requires that the up and down quarks of the third color remain un- 
gapped; this requirement must, therefore, be modified once. these quarks pair 
to form a J = 1 condensate, breaking rotational invarianceH 

3.2 Two+One Flavors 

Nature chooses two light quarks and one middle-weight strange quark. If we 
imagine beginning with the CFL phase and increasing rris, how do we get to 
the 2SC phase? This question has been answered in Refs. 13,17. A nonzero 
rUs weakens those condensates which involve pairing between light and strange 
quarks. The CFL phase requires nonzero (us) and (ds) condensates; because 
these condensates pair quarks with differing Fermi momenta they can only exist 
if the resulting gaps (call them A^^ and A^is) are larger than of order m^/2^, 
the difference between the u and s Fermi momenta in the absence of pairing. 
This means that as a function of increasing TOs at fixed fi (or decmasing /i 
at fixed m-s) there must be a first order unlocking phase transitionfijO The 
argument can be phrased thus: the 2SC and CFL phases must be separated by 
a phase transition, because chiral symmetry is broken in the CFL phase but not 
in the 2SC phase; suppose this transition were second order; this would require 
and A^^ to be infinitesimally small but nonzero just on the CFL side of 
the transition; however, these gaps must be greater than of order rn^/2/x; a 
second order phase transition is therefore a logical impossibility, either in mean 
field theory or beyond; the transition must therefore be first order. Note that 
the rris that appears in these estimates is a density dependent effective strange 
quark mass, somewhat greater than the current quark mass. 

Putting in reasonable numbers for quark matter which may arise in com- 
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pact Stars, for = 200 - 300 MeV and /i = 400 - 500 MeV we find that the 
CFL phase is obtained if the interactions are strong enough to generate a gap 
A which is larger than about 40 — 110 MeV, while the 2SC phase is obtained if 
A is smaller. As we shall see in the next section, A 40 — 110 MeV is within 
the range of current estimates and present calculational methods are therefore 
not precise enough to determine whether quark matter with these parameters 
is in the CFL or 2SC phases. At asymptotically high densities, however, the 
CFL phase is necessarily favored. 

3.S One Flavor 

The remaining variant of color superconductivity which may arise in nature 
is that involving pairing between quarks of a single flavor. If nis is too large 
for CFL pairing (i.e. to^ greater than of order y/2fiA) but is nevertheless 
less than n itself, the 2SC phase will include a nonzero density of strange 
quarks. Although, by hypothesis, these strange quarks do not pair with the 
light quarks, they can pair with each other, forming a (ss) condensate with 
angular momentum J = 1. This pairing has been analyzed in Ref. 47. (See 
also Ref. 48.) Interestingly, even though the Cooper pairs have J = 1, the 
condensate does not break rotational invariancecZI The Cooper pairs are in 
the color 3 channel, and the condensate locks color and spatial rotation: it 
leaves unbroken a global symmetry assocated with simulataneous color and 
spatial rotation. The resulting gap is much smaller than the J 2SC and 
CFL gaps. It has been estimated to be of order hundreds of keVpZl although 
applying results of Ref. 49 suggests a somewhat smaller gap, around 10 keV. 

3.4 Four or More Flavors 

We end this section with brief mention of four variants which are unphysical, 
but nevertheless instructive: QCD with more than three light flavors, QCD 
with two colors, QCD with many colors, and QCD with large isospin density 
and zero baryon density. 

Dense quark matter in QCD with more than three flavors was studied 
in Ref. 10. The main result is that the color-flavor locking phenomenon per- 
sists: Condensates form which lock color rotations to flavor rotations, and 
the SU{Nf)L X SU{Nf)ji group is broken down to a vector subgroup. Unhke 
with Nf =^ 3, however, the unbroken group is not the full SU{Nf)L^fi^ which 
is unbroken in the vacuum. In the case of Nf — 4, for example, one finds 
S'[/(4) L X -S£{(4)i? ^ 0(4) L+R while in the case of Nf = 5, SU{5) lxSU{5)r^ 
5'C/(2)L+ijllj For Nf = 4, 5 as for Nf — 3, chiral symmetry is broken in dense 
quark matter. However, because the unbroken vector groups are smaller than 
SU{Nf)v , there must be a phase transition between hadronic matter and dense 
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quark matter in these theoriesli3 

If Nf is a multiple of three, the order parameter takes the form of multiple 
copies of-jthe Nf = 3 order parameter, each locking a block of three flavors 
to colorty All quarks are gapped in this phase, as in the Nf ~ 3 CFL phase. 
For Nf = 6, the resulting symmetcy breaking pattern is SU{6)l x SU{6)r 
SUi3)L+R X Uil)L+R X f/(l)L-fllia The unbroken SU{3)l+r is a simultaneous 
rotation of both three flavor blocks for L and R and a global color rotation. 
Note that the unbroken [/(l)'s are subgroups of the original SU{6) groups: 
they correspond to vector and axial flavor rotations which rotate one three 
flavor block relative to the other. Note that for Nf — 6, unlike for Nf — 3, 4, 5, 
chiral symmetry is not completely broken at high density: an axial U{1) sub- 
group remains unbroken. As the primary condensate we have just described 
leaves no quarks ungapped, there is no reason to expect the formation of any 
subdominant condensate which could break the unbroken chiral symmetry. 
Both because of this unbroken chiral symmetry and because the unbroken vec- 
tor symmetry differs from that of the vacuum, there must be a phase transition 
between hadronic matter and dense quark matter in QCD with Nf = 6113 

3.5 Two Colors 

The simplest case of all to analyze is QCD with two colors and two flavors. 
The condensate is antisymmetric in color and flavor, and is therefore a singlet 
in both color and flavor. Because it is a singlet in color, dense quark matter in 
this theory is not a color superconductor. Although the condensate is a singlet 
under the ordinary SU{2)l x SU{2)ji flavor group, it nevertheless does break 
symmetries because the symmetry of the vacuum in QCD with Nf — Nc = 2 
is enhanced to SU{4:). One reason why Nc = 2 QCD is interesting to study 
at nonzero density is Aat it provides an example where quark pairing can be 
studied on the latticeE3 The Nc = 2 case has also been studied analytically in 
Refs. 6,51; pairing in this theory is simpler to analyze because quark Cooper 
pairs are color singlets. We refer the reader to these references for details. 

3.6 Many Colors 

The Nc ^ oo limit of QCD is often one in which hard problems become 
tractable. However, the grauiid state of TVc = oo QCD is a chiral density wave, 
not a color superconductorE3 At asymptotically hj^h-jdensities, color supercon- 
ductivity persists up to A'c's of order thousands Eja before being supplanted 
by the phase described in Ref. 52. At any finite N^ color superconductivity 
occurs at arbitrarily weak coupling whereas the chiral density wave does not. 
For Nc = -i, color superconductivity is still favored over the chiral density 
wave (although not by much) even if the interaction is so strong that the color 
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3.7 QCD at Large Isospin Density 

The phase of A^c = 3 QCD with nonzero isospin densifai. (/i/ ^ 0) and zero 
baryon density (/i = 0) can be simulated on the latticeE^ The sign problems 
that plague simulations at ^ do not arise for /z/ ^ 0. Although not physi- 
cally realizable, physics with fij ^ and /i = is very interesting to consider 
because phenomena arise which are similar to those occurring at large and, in 
this context, these phenomena are accessible to numerical "experiments" . Such 
lattice simulations can be used to test calculational methods which have also 
been applied at large ^, where lattice simulation is unavailable. At low isospin 
density, this theory describes a dilute gas of Bose-condensed pions. Large 
fii physics features large Fermi surfaces for down quarks and anti-up quarks. 
Cooper pairing of down and anti-up quarks, and a gap whose g-depeadence 
is as in ( |108D , albeit with a different coefficient of l/g in the exponentcS This 
condensate has the same quantum numbers as the pion condensate expected at 
much lower /x/ , which means that a hypothesis of continuity between hadronic 
— in this case pionic — and quark matter as a function of /x/. Both the dilute 
pion gas limit and the asymptotically large /i/ limit can be treated analyti- 
cally; the possibility of continuity between these two limits can be tested on the 
latticeE3 The transition from a weak coupling superconductor with condensed 
Cooper pairs to a gas of tightly bound bosons which form a Bose condensate 
can be studied in a completely controlled fashion. 



4 Calculational Methods 

In this Section we review the basic theoretical methods used to analyze color 
superconductivity. To keep the discussion within appropriate bounds, we shall 
concentrate on the methods used to characterize the ground state. The deriva- 
tion of effective theories for the low-energy dynamics, electromagnetic response, 
and transport properties is of course predicated on identification of the ground 
state, but presents many additional points of interest as we have already seen 
in Section 2. In that Section, we deferred all discussion of methods by which 
the gap A, which characterizes the ground state, is calculated. Much effort 

prasS|M&3fflSS^^fiiffl§6lfillfflHliffand^n this^sec- 
tion we face up to the challenge of describing what has been learned. 

It would be ideal if the calculation of the gap were within the scope of 
lattice gauge theory as is, for example, the calculation of the critical tempera- 
ture on the vertical axis of the phase diagram. Unfortunately, lattice methods 
relying on importance sampling have to this point been rendered exponen- 



The Condensed Matter Physics of QCD 



27 



tially impractical at nonzero baryon density by the complex action at nonzero 
/J,. Various lattice methods can be applied for /i 7^ as long as T / /i is large 
enough;Ld so far, though, none have proved applicable at temperatures which 
are low enough that color superconductivity occurs. As we saw in the previ- 
ous section, lattice simulations are possible in two-color QCD and in QCD at 
large isospin density. Finally, sophisticated algorithms have recently allowed 
theories which are simpler than QCD but which have as severe a fermioiL^ign 
problem as that in QCD at nonzero chemical potential to be siniulatedllj All 
of this bodes well for the future. 

To date, in the absence of suitable lattice methods, quantitative analy- 
ses of color superconductivity have followed two distinct strategies. The first 
approach is utilitarian and semi-phenomenological, emphasizing the use of sim- 
plified models. This will occupy us in Sections 4.1-4.3. The overarching theme 
here is to define models which incorporate the salient physical effects, yet are 
tractable using known mathematical techniques of quantum many-body the- 
ory. Free parameters within a model of choice are chosen to give reasonable 
vacuum physics. Examples include analyses in which the interaction between 
quarks is replaced simply by four-femiion interactions with the quantum num- 
bers of the inatanton interaction BQ^a or of one-gluon exchangePii3 random 
matrix modeda _and more sophisticated analyses done using the instanton 
liquid modelEJuH Renormalization group methods have also beeii|_used to 
explore the space of all possible effective four-fermion interactions EZIl3 These 
methods yield results which are in qualitative agreement: the favored conden- 
sates are as described in Sections 2 and 3; the gaps range between several tens 
of MeV up to of order 100 MeV; the associated critical temperatures (above 
which the diquark condensates vanish) can be as large as about Tc ~ 50 MeV. 
This agreement between different models reflects the fact that what matters 
most is simply the strength of the attraction between quarks in the color 3 
channel, and by fixing the parameters of the model interaction to fit, say, the 
magnitude of the vacuum chiral condensate, one ends up with attractions of 
similar strengths in different models. 

The second, more ambitious approach is fully microscopic. Such an ap- 
proach is feasible, for high-density QCD, due to asymptotic freedom. Several 
important results have been obtained from th£_microscopic approach, perhaps 
most notably the asymptotic form of the gapEJ Very significant challenges re- 
main, however. It is not really known, for example, how to calculate corrections 
to the leading term in a systematic way. We review the microscopic approach 
in Sections 4.4-4.5. 

These approaches have complementary virtues - simplicity versus rigor, 
openness to phenomenological input versus quantitative predictive power at 
asymptotically high density. Fortunately, they broadly agree in their main 
conclusions as to the patterns of symmetry breaking and the magnitude of the 
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gap at accessible densities of order fx = 400 — 500 MeV. 
4-1 Renormalization Toward the Fermi Surface 

The adequacy of weak coupling for describing QCD at high density is by no 
means obvious. Our strategy will be to adopt it as a working hypothesis, 
bring out its consequences, and see whether we find a consistent picture. In 
doing this, we must consider how to work from the fundamental interactions 
in the perturbative, no-particle state ("Fock vacuum"), which are assumed to 
be simple, to the effective interactions in the dense medium. 

To begin, let us focus on the quark degrees of freedom. In line with our 
announced strategy, we begin by approximating the ground state with filled 
Fermi spheres for all the quarks, with Fermi momentum pp. The low-energy 
excitations then include states where some modes just below the nominal Fermi 
surfaces are vacant (hole modes) and some modes just above are occupied (par- 
ticle modes). There is a continuum of such states. Moreover, states containing 
pairs of particle or hole modes with three-momenta (p, —p) , for various values 
of the direction of p but a common \p\, are all degenerate. We are therefore 
faced with a system whose excitations have a large density of states (propor- 
tional to the area of the Fermi surface) which can be excited at arbitrarily 
small free energy cost (for pairs whose \p\ is arbitrarily close to pp-) In the 
absence of interaction, there would be nothing more to say. The effect of inter- 
actions, however, is to allow all pairs (p, —p) with a common \p\ to scatter into 
one another, consistent with momentum conservation. Thus, we are perturb- 
ing a system with a continuum of excitations with arbitrarily low energies. In 
this situation, with large numbers of nearly degenerate states, straightforward 
perturbation theory can fail, even for weak coupling. 

The Wilsonian renormalization group is often an appropriate tool for ana- 
lyzing problems of this sort. Following this approach, one attempts to map the 
original problem onto a problem with fewer degrees of freedom, by integrating 
out the effect of the higher-energy (or, in a relativistic theory, more virtual) 
modes. Then one finds a new formulation of the problem, in a smaller space, 
with new couplings. In favorable cases the reformulated problem is simpler 
than the original, and one can turn around and solve it. 

This use of the renormalization group is in the same spirit as that ap- 
propriate for the analysis of long wavelength physics at a second order phase 
transition, or the long wavelength physics of the QCD vacuum (a.k.a chi- 
ral perturbation theory.) In traditional perturbative QCD one nms the logic 
backwards. The fundamental short-distance theory, with nothing integrated 
out, is simple and weakly coupled. When one integrates out highly virtual 
modes, one finds that QCD becomes more strongly coupled. Thus, renormal- 
ization group methods are used in two distinct ways in vacuum QCD: first. 
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they inform us how the fundamentally simple theory of quarks and gluons 
comes to look complicated at low energy and, second, once rephrased in terms 
of the correct effective degrees of freedom relevant at low energies, they allow 
us to make predictions for low energy physics. Renormalization group methods 
are used in two analogous ways at high density. In this Section we shall use 
them to inform us at what free energy scale the theory based on the original 
quark degrees of freedom breaks down, signalling the formation of a gap and 
the need to switch to new degrees of freedom. These new degrees of freedom 
include the gapped quasiparticles, which can then be integrated out. If the 
gapped phase features massless or light Nambu-Goldstone boson excitations of 
the condensate then, as we have seen in Section 2, we can use renormalization 
group methods to make predictions for the physics of these effective degrees of 
freedom at momenta well below the gap. 

In both chiral perturbation theory of the QCD vacuum and the Wilsonian 
effective theory appropriate at high density, one integrates out modes with 
large free energy, keeping only those with smaller and smaller free energy. In 
the vacuum, free energy is just energy and this corresponds to keeping modes 
with momenta closer and closer to zero. At high density, however, the modes 
with small free energy are those with momenta near the Fermi surface, and it 
is on these modes that the renormalization group focuses. The application rif 
Wilsonian methijds to physics at a Fermi surface was pioneered by ShankaO 
and Polchinskio 

Asymptotic freedom plays a foundational role in the analysis of QCD at 
asymptotically high density. It allows us to make it plausible that the funda- 
mental couplings among the active modes, near the Fermi surface, are gener- 
ically weak. Indeed, these modes involve large energy and momenta, and so 
scattering events which displace them by a fixed finite angle involve large mo- 
mentum transfer, whereas scattering events with small angle should not much 
alter the physical properties of the state. This intuitive argument will be re- 
fined, and justified, in Section 4.4. 

Granting that the fundamental coupling is weak, let us consider the con- 
sequences of integrating out modes whose free energy is between e and 5e on 
the effective interactions among the remaining modes, which are those within 
a band in momentum space centered on the Fermi surface which have free en- 
ergy less than 5e. We will obtain renormalized couplings among the remaining 
modes, due to one-loop diagrams in which external legs are within 5e of the 
Fermi surface while loop momenta are between 5e and e. This approach was 
first applied to high density QCD in Ref. 57 and was generalized in Ref. 58. 

The first result of this analysis is that, in general, four-fermion, six-fermion 
and highetorder interactions are all suppressed as we approach the Fermi 
surfacellj'Lj This fixed point corresponds to Landau Fermi liquid theory. The 
only, and crucial, exceptions are those four-fermion operators that involve in- 
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volve particles or holes with equal and opposite three momenta. These allow 
pairs of particles on the Fermi surface to scatter, as mentioned above. For 
couplings G,j of this kind we find 

d]^5 " " ""^^ ■ ^^^^ 
Here is a catch-all label that might include color, flavor, angular momentum. 
The effective couplings will depend on all these indices. In general we will have 
a matrix equation, with couplings between the different channels (different 
values of 77). But to bring out the main point, it is enough to consider the 
simplest possible case, with one channel, as we have done in (|l5|). In this 
single channel calculation, a simple calculation shows that n is positive. (Note 
that we have chosen a notation in which attractive interactions correspond to 
G > 0.) This evolution equation is quite simple to integrate, and we find 

Thus we see that the qualitative behavior of G,j{S) depends on the sign of 
Gj^{l). If this is negative, then as (5 ^ one finds that G,^{6) from below. 
Repulsive interactions are thus irrelevant. On the other hand, if G,,(l) > 
(indicating an attractive coupling) then Gr,{S) grows, and formally it diverges 
when 

^ = ^'^P ( • (1^) 

V K^G.,j{l)J 

Of course once the effective coupling becomes large the working equation (^, 
which neglects higher-order corrections, is no longer trustworthy. All we may 
legitimately infer from (|lj) is that an attractive interaction in any channel, 
however weak, will, by renormalization toward the Fermi surface, induce a 
strong effective coupling between particles with equal and opposite momenta, 
on opposite sides of the Fermi surface. 

What we have done here is simply to rephrase Cooper's discovery of the 
pairing instability in modern language, suitable for generalization. Now let us 
apply it to four-fermion interactions of the type we expect to occur in QCD. 

We will restrict ourselves to massless QCD, a spherical Fermi surface, and 
local operators invariant under the appropriate chiral symmetry. The cases of 
three and two ffavors are rather different, and we will analyze them separately 
in turn. 

Following Ref. 57, we take the basic four-fermion interactions in Nf = 3 
QCD to be proportional to the operators 

OIl = {i'L7oM\ OIr = {i^LlO^LK^BHoi^R) , (18) 
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Each of these operators comes in two color structures, for example color sym- 
metric and color anti-symmetric 

{SabScd ± 5adSbc) ■ (19) 

Nothing essentially new emerges upon considering superficially different isospin 
structures, or different Dirac matrices. All such structures can be reduced to 
linear combinations of the basic ones (^8|), or their parity conjugates, by Fierz 
rearrangements. In total, then, we need to consider eight operators. 

As discussed above, these operators are renormalized by quark-quark scat- 
tering in the vicinity of the Fermi surface. This means that both incoming and 
outgoing quarks have momenta pi,P2 — ±P and P3,P4 — ±g with |g| — Pf- 
We set the_j3xternal frequency to be zero. A graph with vertices Fi and F2 
then givesEJ 



G1G2I (Fi),,,(Fi)fc,fc 



(F2),,'(F2),^ (20) 



with / — g^//^ log((5). We will denote the density of states on the Fermi 
surface by = /i^/(27r^) and the logarithm of the scale as t ~ log{6). The 
renormalization group does not mix LL and LR operators, nor different color 
structures. This means that the evolution equations contain at most 2x2 
blocks. A simple calculation now yield£Zl 



dt 3 
d{G(i^ - 3Gf ^) _ ^^^^^^ 
dt 

d{G^^ + iO^") 
Jt 

d{GQ^ - Gf^) _ 2A^ ^LB\2 



(G^^ + Gt^r , (21) 

Ar(G^^ - 3Gp)2 , (22) 

, (23) 

(Go^^ - Gf^f . (24) 



dt 3 

In this basis the evolution equations are already diagonal. The coupling Gi = 
Gl^^ + Gf^ evolves as 

^'^^^ = l + (A^/3)Gi(0)t' ^^^^ 

with analogous results for the other operators. Note that the evolution starts 
at t = and moves towards the Fermi surface as t — s- — 00. If the coupling is 
attractive at the matching scale, Gi(0) > 0, it will grow during the evolution, 
and reach a Landau pole at t^ = 3/{NGi{0)). The location of the pole is 
controlled by the initial value of the coupling and the coefficient in the evolution 
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equation. If the initial coupling is negative, the coupling decreases during the 
evolution. The second operator in ( pT| ) has the largest coefficient and will reach 
the Landau pole first, unless the initial value is very small or negative. In that 
case, subdominant operators can come to cause the leading instability. 

The form oLthe operators that diagonalize the evolution equations is read- 
ily understood^ First, order the operators according to the size of the coeffi- 



cient in the evolution equations 

Odom = {^LlO^l^Lf - {i^Ll^L? , (26) 

Osub,i = {4'l1o'4^l){'4'r1o^r) - ^{^-'l1^l){4'r1^r) , (27) 

Osub,2 = {i^Lloi^L? + \{i'Ll^L? , (28) 

Omar = {i'Lloi^L){i'Rl()-4^R) + {■4'L1'iPl){-4'R1^r) ■ (29) 

Upon Fierz rearrangement we find 

Odom = 2{i^LCi^L){i^LCi^L). (30) 

Osub,l = \{i^LCli^R){^RCl^L) + . ■ ■ , (31) 

Osub,2 = ^{'iPLCt^L){^LCt^L). (32) 

Omar = ^ (V'lC'TO V'i? ) (V'i? V'L ) + • • ■ • (33) 



This demonstrates that the linear combinations in ( p6[ - |29| ) correspond to simple 
structures in the quark-quark channel. (It also means that it might have been 
more natural to perform the whole calculation directly in a basis of diquark 
operators.) 

The full structure of the {LR) operators is 

Osub,l,Omar = (V'C'7Ts^ A "0) ('/'CT'^S, aV") + (V'C'775TA,S'/') ("0^775X^,5 V^) , 

where ts^a are symmetric, respectively anti-symmetric, flavor generators. Note 
that because these two structures have different flavor symmetry, the flavor 
structure cannot be factored out. 

The dominant operator corresponds to the scalar diquark channel, while 
the subdominant operators contain vector diquarks. Note that one cannot 
decide which color channel is preferred from the evolution equation alone. To 
decide that question, we can appeal to the fact that "reasonable" interactions, 
including specifically one gluon exchange, will be attractive in the color anti- 
symmetric but repulsive in the color symmetric channel. Indeed, it is the color 
anti-symmetric configuration that minimizes the total color flux (and hence 
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field energy) emanating from the quark pair. If the color wave function is 
anti-symmetric, the dominant operator fixes the isospin wave function to be 
anti-symmetric as well. 

The form of the dominant operator indicates the existence of potential 
instabilities, but does not itself indicate how they are resolved. The renormal- 
ization group analysis tells us the scale at which the dominant couplings get 
strong but, strictly speaking, tells us nothing more. In Sections 4.2 and 4.3, 
we shall see how to learn more, and in particular to confirm that the resolution 
of the instability is that a gap forms in the dominant attractive channel. Be- 
cause the renormalization group methods do not allow us to directly describe 
the formation of the gap (although they indicate its magnitude) they cannot 
be used to confirm that color-flavor locking is favored. For this, a variational 
calculation is required as we shall describe in Section 4.2. 

The dominant four-fermion operator in three-flavor QCD does not distin- 
guish between scalar and pseudoscalar diquarks. Indeed, for Nf > 3 all four- 
quark operators consistent with chiral symmetry exhibit an accidental axial 
baryon symmetry, under which scalar and pseudoscalar diquarks pSJi^e-equiva- 
lent. For Nf = 3 this degeneracy is lifted by six-fermion operatorsPtl'td which 
are irrelevant operators in the sense of the renormalization group. In the mi- 
croscopic theory, these operators will be induced by instantons. Their formal 
irrelevance does not imply that they are negligible physically, especially since 
they are the leading terms which break the residual axial baryon number sym- 
metry. 

Anticipating that the instability indicated by renormalization toward the 
Fermi surface is stabilized by the formation of a gap, we conclude that upon 
further renormalization the couplings cease to run. We therefore expect that, 
instead of running to zero, the instanton coupling remains at some finite 
value. Instantons have important physical effects in the condensed phase, even 



antiquark pairs to condense, by inducing the "normal" chirality-violating order 
parameter (familiar at zero density) from the diquark-antidiquark condensate 
(which arose, in turn, as a secondary consequence of the primary diquatk con- 
densation). The resulting quark- antiquark condensate is very smallla't2l The 
instanton interaction also lifts the degeneracy between the scalar and pseu- 
doscalar diquark condensates, favoring the scalar. 

We turn now to QCD with two flavors. This requires us to take into 
account additional operators. At first hearing it might seem odd that with 
fewer flavors we encounter more possible interactions. It occurs because for 
Nf = 2, but not for larger values, two quarks of the same chirality can form 
a chiral 317(2) x SU(2) singlet. Related to this, for Nf ^ 2 we have U{1)a 
violating /our- fermion operators. In the microscopic theory, these operators 
will be induced by instantons. 




Specifically, they cause quark- 
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The new operators are 

Os = deii^Ril^L), Ot = det(V5flSVL) ■ (34) 

Both operators are determinants in flavor space. For quark-quark scattering, 
this impUes that the two quarks have to have different flavors. The fact that 
the flavor structure is fixed implies that the color structure is fixed, too. For 
a given [qq) spin, only one of the two color structures contributes. Finally, 
both quarks have to have the same chirality, and the chirality is flipped by the 
interaction. 

These considerations determine the structure of the evolution equations. 
Two left handed quarks can interact via one of the instanton operators, become 
right handed, and then rescatter through an anti-instanton, or through one of 
the U{1)a symmetric RR operators. The result will be a renormalization of 
the LL vertex in the first case, and a renormalization of the instanton in the 
second. Note that the flavor structure will always remain a determinant. Even 
though instantons generate all possible Dirac structures in (^8|), the color- flavor 
structure is more restricted. 

Evidently, instantons do not affect the evolution of the LR couplings at 
all. The evolution equations of the LlL couplings are modifled to become 
(henceforth we drop the subscript LL) : E3 





N 


dt 




dOi 


N 


dt 




dKs 


N 


dt 




dKr 


N 


dt 


Y 



2GaGi - hGj -Kl + 2KsKt - 5K^ } , (35) 




- ^ 2 (-Go + G,) Ks + 2 {Go- 5G,;) Kt 



^ (Go - 5G,) Ks+'^ (-5Go + 13G,) Krj • (38) 
These equations can be uncoupled in the form 

^ = -N{Gl+Kl), (41) 



dK. 



2 



dt 



-2N G2K2 , (42) 



The Condensed Matter Physics of QCD 



35 



where Gi = Go + G„ Ki = Ks + Kt and G2 = Go - 3Gj, K2 ^ Kg - iKr- 
The equations for G, i^T uncouple even further. We have 



d{G2 + K2) 
dt 

d{G2 - K2) 
dt 



-N (G2 + K2) 
-N (G2 - K2) 



as well as the analogous equation for Gi , Ki. 

The differential equations are now easily solved^ leading to 



G2{t) = - 



1 



1 



2 [a + Nt b + Nt 



K2{t) = - 



1 



1 



2 \a + Nt b + Nt 



(43) 
(44) 

(45) 
(46) 



together with the analogous results for Gi, Ki. Here, a,b = (G2(0)±_ft'2(0))~"'^. 
We see that G2 and K2 will grow and eventually reach a Landau pole if either 
a or is positive. The location of the pole is determined by the smaller of the 
values, tc = —a/N or tc = —b/N. The same is true for Gi and Ki, but the 
couplings evolve more slowly, and the Landau pole is reached later. 

At this level a number of qualitatively different scenarios_are possible, 
depending on the sign and relative magnitude of G(0) and if (0)L3 (Henceforth 
we drop all subscripts.) If G(0) and K{0) are both positive then they will both 
grow, and the location of the nearest Landau pole is determined by G{0)+K{0). 
The asymptotic ratio of the two couplings is 1. If G(0) and K{Q) are both 
negative, and the magnitude of G(0) is bigger than the magnitude of K(Q), 
then the evolution drives both couplings to zero. These are the standard cases. 
Attraction leads to an instability, and repulsive forces are suppressed. 

More interesting cases arise when the sign of the two couplings is-differ- 
ent. The case G(0), A'(0) < and \K{0)\ > |G(0)| is especially weirdN Both 
G(0), A'(0) are repulsive, but the evolution drives G(0) to positive values. Both 
couplings reach a Landau pole, and near the pole their asymptotic ratio ap- 
proaches minus one. Similarly, we can have a negative G(0) and positive KLQ^ 
with K{0) > |G(0)|. Again, the evolution will drive G(0) to positive valuesE3 

The dominant and sub-dominant instanton operators are 



Gdom 
O suh 



det 
/ 

det 
/ 



{Ml? + 7;i^R^i^L)^ 



Upon Fierz rearrangement, we find 

Odom = 2(V'LGT2l/'L)(V'i?CT2?/5fl) 



(47) 
(48) 

(49) 
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Osub = -[^LCT2ti^L){i^RCT2ti^R) , (50) 

corresponding to scalar and tensor diquarks. Both operators are flavor singlet. 
Overall symmetry then fixes the color wave functions, anti-symmetric 3 for 
the scalar, and symmetric 6 for the tensor. The dominant pairing induced by 
instantons is in the scalar diquark channel, the only other attractive channel 
is the tensor. All this neatly confirms the scenario discussed in Refs. 5,6 and 
Section 3.1. 

Just as we found for Nf > 3, there is an appealing heuristic understand- 
ing for the amazingly simple behavior of the evolution equations, obtained 
by focussing on the diquark channels. Instantons distinguish between scalar 
diquarks with positive and negative parity. G + K corresponds to the posi- 
tive parity operator (ipCj^ip) and G — to the negative parity {ipCtlj). The 
asymptotic approach of G/K —f 1, then corresponds to the fact that scalar 
diquark condensation is favored over pseudoscalar diquark condensation. This 
is always the case if K{0) > 0. 

We can also understand the strange case G{0),K{0) < and |i^(0)| > 
|G(0)|. In this case the interaction for scalar diquarks is repulsive, but the in- 
teraction in the pseudoscalar channel is attractive and leads to an instability. 
Note that this can only happen if we have the "wrong" sign of the instanton 
interaction, i.e. for = tt. Similarly, we can understand why the asymp- 
totic ratio of the molecular (instanton-anti-instanton) and direct instanton 
couplings approaches G/K = ±1. Instantons induce a repulsive interaction 
for pseudoscalar diquarks. During the evolution, this coupling will be sup- 
pressed, whereas the attractive scalar interaction grows. But this means that 
in the pseudoscalar channel, the repulsive (instanton) and attractive (molec- 
ular) forces have to cancel in the asymptotic limit, so the effective couplings 
become equal. 

To match all the instanton coupling constants to microscopic QCD would 
require control of the instanton density, the relevant value of as, and the 
screening mechanism, and does not seem practicable. Some simple qualitative 
conclusions may be inferred, however. From the form of the instanton verte.Y we. 
can fix the ratio of the two instanton-like couplings, Kt/Ks = l/(2iVc— 1)ESC3 
This emphasizes that the tensor channel coupling is expected to be small. 

Also, both one gluon exchange and ideas based on instantons (e.g., the 
instanton liquid model) yield Gq ^ > 0, —Gf^ > 0. Thus the favored scenario 
is that both instanton and U(1)a symmetric couplings fiow at the same rate. 
In any case, the leading instability occurs in a scalar diquark channel. 

The simplified analysis presented here is incomplete, in that we have re- 
tained only instantaneous s-wave interactions. Realistic interactions, such 
as those generated by gluon exchange or instantons, are momentum- and 
frequency-dependent. Nevertheless the simplified analysis is valuable, in that 
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it indicates the internal quantum numbers of the channels likely to be most 
unstable. Moreover, other things being equal, s-wave instabilities seem likely 
to dominate, since they add constructively over the fermi surface. (This intu- 
ition is borne out by variational calculations calculations of gaps and pairing 
energies associated withjJiigher angular momentum condensates, which turn 
out to be very smallBoa) 

4.2 Model Hamiltonian by Variational Methods 

The renormalization group determines the running of couplings from the match- 
ing point downwards, towards the Fermi surface. The catastrophic growth of 
an attractive interaction we have encountered indicates an instability of the 
system, indicates in what channel(s) the instability is most pronounced, and 
indicates the energy scale at which the instability develops. At the scale at 
which the instability sets in, however, the perturbative renormalization proce- 
dure bceaks down. The great achievement of Bardeen, Cooper, and Schrieffer 
(BCSpJ was to demonstrate how Fermi-surface instabilities of this type are 
resolved by pairing, to produce a superconducting ground state. This ground 
state is characterized by a gap A, whose magnitude does indeed turn out to be 
of the same order as the energy scale indicated by the renormalization group 
analysis. The instability arose in a renormalization group analysis done in 
terms of fermionic excitations (about the naive ground state with no pairing) 
which have zero free energy at the Fermi surface. In the gapped state, which 
the BCS analysis determines is formed, the resolution of the instability is sim- 
ply that the correct fermionic degrees of freedom are quasiparticle excitations 
which have a gap in their spectrum. In this section we outline the BCS ar- 
gument that the resolution of the instability is a paired state, and give the 
associated calculation of A. 

As an aside, note that once A is known and the ground state has been 
characterized, if (as in the CFL phase) there are bosonic excitations of this 
ground state with energy less than A, a renormalization group analysis framed 
in terms of these new degrees of freedom can be employed in their description. 
The results of such analyses were presented in Section 2.3. 

The original BCS theory was founded on a variational calculation using 
a model Hamiltonian. Their model Hamiltonian retained only an idealized 
form of the effective attraction near the Fermi surface, that encouraged con- 
densation in the ajvave, spin-singlet channel. The pioneering work of Nambu 
and Jona-LasinioEj who first applied the ideas behind BCS theory to particle- 
anti-particle pairing in quantum field theories describing elementary particle 
interactions, followed a similar approach. We shall follow the literature in re- 
ferring to models of this class as NJL models, but shall apply them a la BCS, 
to quark-quark pairing. 
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In QCD there are many more fermion species (colors and flavors) than in 
conventional superconductors, and the form of the condensate is a major focus 
of interest. Renormalization group calculations, as presented in Section 4.1, 
provide important guidance. The dominant operator indicates the existence 
of potential instabilities that are amplified by renormalization, but does not 
indicate how they are resolved. To decide the most favorable form of pairing, 
we must compare the energies associated with different possibilities. For this 
qualitative purpose, and of course to estimate the quantitative relationship 
between the size of the pairing gap and the microscopic parameters, we still 
require detailed dynamical calculations. NJL models provide an appropriate, 
tractable starting pointurEl One chooses a simple four-quark interaction, which 
provides attraction in the channel which the renormalization group analysis 
indicates is dominant. One first normalizes the parameters in such a toy model 
to give reasonable vacuum physics, and then uses the model to estimate the 
gap. Here we shall exemplify the use of such models by using them to estimate 
the size of the gaps in the 2SC and CFL phases and to argue that the CFL 
condensate is indeed the favored pattern of condensation in three-flavor QCD. 

We use an NJL model with free energy 

n^fd^xi,{x){f-ti-foMx)+Hj, (51) 



where the interaction Hamiltonian 

Hi^\g j SxT ii^ixh^T^^ix)) (i^ixh'^T^^ix)) (52) 

describes a four-fermion interaction with the color, flavor, and spinor structure 
of single-gluon exchange. 

In microscopic QCD the interactions become weak at high momentum. 
Here this feature is caricatured with a form factor J-. In detail, when we 
expand Hi in momentum modes we are instructed to include a form factor 
F{p) on each leg of the interaction vertex. Examples which have b£en_used in 
the literature include power-law and smoothed-step profiles for F:EnZl'tj 

or F(p)=(^l + cxp[^]) . (53) 

One of the important results found in these investigations is that if, for any 
given F{p), the coupling constant G is chosen to yield a reasonable vacuum 
chiral condensate (say with a constituent quark mass of 400 MeV) then the 
resulting superconducting gap is reasonably insensitive to the choice of F{p). 
Varying the scale A and the shape of the form factor while tuning the coupling 
to keep some physical property of the vacuum fixed affects the gap very little. 
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This robustness indicates that the magnitude of the gap is largely determined 
just by the strength of the interaction (normalized with respect to a physical 
observable) and does not depend sensitively on other details. For this reason, 
and following Ref. 13,49 we make the simplifying choice of replacing the smooth 
form factor F{p) by a step function: Fijj) = 1 for p < K and F{p) = for 
p > A. 

We begin with two flavors and follow the analysis presented in Refs. 5,49. 
We seek a condensate of the 2SC form by making the BCS ansatz 

I*) = 44J0), 

4 = n H^i(P) + ^"'''e'^'^''^«i^^i(P)«LjP)«U(-P)) ' (54) 

p,Q,/5 

^Ij, = as above, L ^ R, 

for the ground state wave function. Here, a and /3 are color indices, u and d 
label flavors explicitly, and is the particle creation operator (for example, 
'^Lda creates a left-handed down quark with color a). The Cooper pairs de- 
scribed by this ansatz are evidently antisymmetric in color; that they are also 
antisymmetric in flavor follows from the anticommutation relations satisfied by 
the creation operators. The 9's and ^'s are the variational parameters of our 
ansatz: they are to be chosen to minimize the free energy of this BCS state. 
Note that we really should have included pairing among antiparticles in our 
ansatz also. However, this doubles the length of the equations and makes little 
difference in the end. We therefore leave the antiparticles out until the end of 
the derivation, and restore them in the gap equation itself. Note also that |0) 
is the no-particle state. At nonzero chemical potential, in the absence of any 
interactions, the ansatz ( |5^ ) describes filled Fermi seas if Ol{p) = 0Rip) = vr/2 
for IpI < fi, and ^l(p) = ^fl(p) = otherwise. Throughout this section, we 
assume /i is the same for all flavors of quarks. The consequences of relaxing 
this assumption are described in Ref. 49 and sketched in Section 6. 
Explicit computation demonstrates that the condensate 

Tl ^ -^me,,ea,p,^^'^{r)CLi,^^{r)\M>) (55) 

is nonvanishing. Here, C — 17^7^ and L — (1 — 75)/2 is the usual left-handed 
projection operator. T]^ can be expressed in terms of the variational parameters 
as: 

Fl - |^sin^?i(p)cos^L(p)e'(«^(P)-^(P» . (56) 
p 

Here V is the spatial volume of the system and the dependence on the az- 
imuthal angle 4> follows from our use of the spinor conventions described in 
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Refs. 5,7,79. The expression for Tr is the same as that in ( p6| ) except that 
0{p) is replaced by tt — (/)(p). In Eq. ( p6| ) and throughout, (l/T^)X]p be- 
comes / cPp/{2tt)^ in an infinite system. Our simphfied choice of form factor 
corresponds to restricting the range of this integral to |p| < A. 

We must now minimize the expectation value of the free energy ( |5l| ) in 
the state |*): 

{^\nm = 4^ (IpI - ^^) sin^ ^^(p) + 4^ (|p| - /i) sin^ ^^(p) + {Hj) (57) 
p p 

where 

CV 

{Hi) = -—{\TL\' + \THn (58) 
in terms of the condensates which are in turn specified in terms of the varia- 



tional parameters by (56). We have neglected the contribution to (fi) made 
by the third color quarks, which do not pair. 

Upon variation with respect to the ^'s, we find that they must be chosen 
to cancel the azimuthal phases 0(p) in (|56|). In this way, we obtain maximum 
coherence in the sums over p, giving the largest possible magnitudes for the 
condensates and gap parameters. We have 

^Lip) ^ ^(p) + , Cfl(p) =7r-0(p) (59) 

where ipL and ipu are arbitrary p-independent angles. These constant phases 
do not affect the free energy — they correspond to the Goldstone bosons for 
the broken left-handed and right-handed baryon number symmetries — and 
are therefore not constrained by the variational procedure. For convenience, 
we set fL^'^B. — O and obtain condensates and gap parameters that are 
purely real. 

The relative phase tpL — fR determines how the condensate transforms 
under a parity transformation. Its value determines whether the condensate 
is scalar, pseudoscalar, or an arbitrary combination of the two. Because single 
gluon exchange cannot change the handedness of a massless quark, the left- 
and right-handed condensates are not coupled in the free energy 17. Our choice 
of interaction Hamiltonian therefore allows an arbitrary choice oi — ^r- A 
global U{1)a transformation changes (^l — 'fR, and indeed this is a symmetry of 
our toy model. If we included [/(l)A-breaking interactions in our Hamiltonian, 
to obtain a more complete description of QCD, we would find that the free 
energy depends on ip^ — ipn, and thus selects a preferred value. For example, 
had we taken Hj to be the two-flavor instanton interaction as in Refs. 5,6, 
the interaction energy would appear as rj^Fj^ -I- F^F^ instead of as in (p8|). 
This would enforce a fixed phaaa^relation (f^ — fR. — 0, favoring the parity 
conserving (ipCj^^) condensateBu 
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We now apply the variational method to determine the angles 0(p) in our 
trial wavefunction, by minimizing the free energy: d{il) /d6{p) ~ 0. Every- 
thing is the same for left and right condensates so we hereafter drop the L and 
R labels. Upon variation with respect to 0(p), we obtain 

tan20(p) = ^— , (60) 

IpI -M 

where /S.l,r — GTl j^. With the 9 angles now expressed in terms of A, we can 
use the wave function \^!) to obtain expressions for the quasiparticle dispersion 
relation: 



i?(p) = V(IP|-/^)' + ^' • (61) 

Here, E is the free energy cost of removing a pair and replacing it with either 
and up quark with momentum p or a down quark with momentum —p. This 
result confirms that A is the gap in the spectrum of fermionic excitations. 
Note, however, that in the 2SC phase the third color ("blue") quarks remain 
gapless in the ansatz within which we are working. As mentioned in Section 
3, it is likely that they form an angular momentum J = 1 condensate with a 
gap which is many orders of magnitude smaller than the gap A for the red and 
green quarks. 

Substituting the expression (^) for the 6 angles into the expression (|5^) for 
the r condensate, and using the relation A = GF, we obtain a self-consistency 
equation for A: 




1 1 



v/(|p| - + A2 V(IpI+m)' + A2 



(62) 



The first term on the right-hand side of the gap equation (|62|) yields a logarith- 
mic divergence at the Fermi surface if A is small. This term is the contribution 
to the gap equation from particles and holes, and the logarithmic divergence is 
the manifestation of the BCS instability. Because the right hand side diverges 
for A — > 0, there must be a solution to the gap equation with A ^ even for an 
arbitrarily small coupling G, as long as G is positive (attractive). The second 
term, with (|p| -l-/i) in the denominator, is the contribution from antiparticles, 
which we have suppressed above and restored here. Note that even though 
the same value of A appears in both the antiparticle and particle/hole contri- 
butions, the antiparticle contribution to the gap equation is small because its 
denominator is everywhere greater than /i. 

We shall rederive the gap equation (^^ diagrammatically in Section 4.3. 
The variational method we have used above has the virtue of providing us in 
addition with an expression for the corresponding wave function itself. 



42 



Handbook of QCD / Volume 3 



In order to use this gap equation to estimate A, we need to fix the coupling 
G. To do this, we evaluate the vacuum chiral gap (or constituent quark mass) 
and set it to M = 400 MeV. In vacuum, M satisfies the gap equation 



which may be derived along very similar lines to the above, working at zero 
chemical potential and postulating a chiral symmetry breaking condensate in 
which quarks and anti-quarks pair in the standard fashion. Note that for 
M ^ 0, the right-hand side of the chiral gap equation has no logarithmic 
divergence; no BCS Fermi surface instability. This means that M 7^ only if 
G is above some threshold. Let us choose A = 800 MeV and fix G such that 
M = 400 MeV. Upon solving the BCS gap equation ( ^ we then find that 
A = 106 MeV when ^ = 400 MeV and A = 145 MeV when = 500 MeV. 

Note that if, instead, we had used the four-fermion interaction with the 
quantum numbers of the instanton vertex, as in Ref. 5, we would have obtained 
precisely the same results. That is, had we made the same choice for the form 
factor and for A, and then chosen the coupling constant for the instanton vertex 
such that the vacuum constituent quark mass is 400 MeV, it turns out that 
the resulting A is exactly that which we have obtained using the single-gluon 
exchange interaction. 

Varying A, or using a smooth form factor instead of a sharp cutoff, yields 
similar resLults for A, ranging from about 50 MeV to somewhat more than 
100 MeVBe^ The authors of Ref. 6 also use the instanton interaction, and also 
normalize it such that M = 400 MeV in vacuum, but they model the fact 
that instanton effects decrease at nonzero density by introducing a density- 
dependent coupling constant. They find gaps ranging from 20 MeV to 90 MeV. 
More sophisticated treatments of the instanton interaction, including form 
factors obtained from suitable Fourier transforms of .instanton profiles, tend 
to yield larger gaps, perhaps as large as 200 MeVOEJ 

The models we are discussing in this subsection can only be used for qual- 
itative guidance, but it is very pleasing to see so many of them agreeing on the 
order of magnitude of the gap. 

The 2SC phase Ima^also been studied at nonzero temperature using the 
instanton interactionBcj The critical temperature above which the condensate 
vanishes (in mean field theory) turns out to be_within a few percent of that 
expected from BCS theory: = 0.57A(r = 0)^1 The standard BCS result is 
valid in the A//x — > limit, and the interaction in QCD is strong enough that 
A//^ ^ 1/4 at the densities we are discussing, so a few percent discrepancy 
is to be expected. Note that based on our estimate of the gap, the critical 
temperature we estimate for the color superconducting phase is of order many 




(63) 
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tens of MeV. This makes it almost certain that this phase is inaccessible in 
heavy ion collisions, which are much hotter. Although Tc is low relative to 
heavy ion collision temperatures, color superconductors are in three senses 
high temperature superconductors. First, ^ 10^^ Kelvin. Second, and more 
seriously, can easily be of order 10% of the Fermi energy. This makes these 
materials high temperature superconductors by any definition, even relative to 
the cuprates. This robustness of the phenomenon of color superconductivity in 
the face of nonzero temperature directly reflects the strength of the attraction 
between quarks in QCD. Third, Tc is much hotter than the temperature in a 
neutron star which is more than a few seconds old, making it clear that if quark 
matter exists in the core of a neutron star, it must be a color superconductor. 

We turn now to QCD with three flavors, and derive a gap equation appro- 
priate for the CFL phase. We use the same single gluon exchange interaction 
Hamiltonian Hj as above. This time, howevec, we present a derivation which 
follows the ideas of Bogoliubov and ValatinE3 We focus from the beginning 
only on the condensate of left-handed quarks. (As above, since we do not in- 
clude any U{1)a breaking interaction the left- and right-handed condensates 
can be rotated one into the other. Introducing instanton effects, done in the 
CFL phase in Refs. 9,10, favors the Lorentz scalar combination.) 

We rewrite the interaction Hamiltonian in terms of Weyl spinors, keeping 
only the left-handed fields: 

Hi^^j cfx^mSl - 5p'l){^^\^%e^'^i^%^'') , (64) 

We have explicitly displayed color {a,(3 . . .), flavor . . .), and spinor (a, d . . .) 
indices, and rewritten the color and flavor generators. 

We now make the mean- field ansatz that in the true ground state \^) at 
a given chemical potential, 

G (65) 
Pl^i^ = i(A8-t-iAi)^^5^-f iAi<5;5i , 

where the numerical factors have been chosen so that Ai and Ag, which param- 
eterize P, will turn out to be the gaps for singlet and octet quark excitations, 
respectively (see below). This ansatz takes the color 3 form (0) if Ai = — 2A8. 
Solutions to the gap equation turn out to be within a few percent of satisfying 
this, meaning that there is only a small admixture of the color 6 condensate. 
In the mean field approximation, the interaction Hamiltonian becomes 

Hi - \jd^xT Q% + c.c. , ^^^^ 

Q% = A85^5^ + i(Al-A8)5^^^ 
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Replacing indices i, (3 with a single color-flavor index p, we can simultaneously 
diagonalize the 9x9 matrices Q and P, and find that they have two eigenvalues, 

Qi = Ai, Q2---Q9 = iAg. 

That is, eight of tlic^ nine quarks in the theory have a gap paramc^tc^r givcm by 
Ag, while the remaining linear combination of the quarks has a gap parameter 
Ai. 

The Hamiltonian can be rewritten in this color-flavor basis in terms of par- 
ticle/antiparticle creation/ annihilation operators ap,bp. (This time, we keep 
the antiparticles throughout the derivation.) We also expand in momentum 
modes (sec Rcfs. 7,79 for our spinor conventions) and now explicitly include 
the form factors F{p), 



H = ^ (ft-/x)aj(k)a,(k)+ ^ (/x-A;)a+(k)a^(k) 

p,k>fj. p,k<fj. 

+ 5^(fc + M)6j(k)6,(k) 
+ iX^F(p)20,e-'*(p)(a,(p)a,(-p) +6t(p)6t(-p)) +C.C. , 

(68) 

where the perturbative ground state, annihilated by Up and bp is the Fermi 
sea, with states up to pi?- = /x occupied. Finally, we change basis to cre- 
ation/annihilation operators y and z for quasiparticles, 



yp{k) = eos(^^;^(k))ap(k)+sin(^^;^(k))exp(^C^(k))at(-k), 

zpik) = cos(^^(k))6p(k)+sin(^^(k))exp(zC^-(k))6t(-k), 

where 

cos(20^(k)) = , e,nk) = <^(k)+7r, 

cos(20^(k)) = , =, C^(k) = -</.(k). 

^(/e + ^)2 + F(fc)4Q2 

(70) 

These values are chosen so that H has the form of a free Hamiltonian for 
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quasiparticles: 



+ ^{k + ^iY+F{kYQl zt(k)z,(k) 



(71) 



Clearly then the ground state j?/;) contains no quasiparticles: 

yp(k)l^) =zp(k)|^) =0. (72) 



The gap equations follow from requiring that the mean field ansatz (65) 
hold in the quasiparticle basis. In other words, we use (^9|) and (|7^) to rewrite 
( |65| ) in terms of quasiparticle creation/annihilation operators, and then evalu- 
ate the expectation value using (^). We get two gap equations: 



where 



i^(A)^4E{ 



Ag + iAi = Gi^(Ai), 
iAi = GA'(A8), 

i^(fc)4A F(fc)4A 



v/(fc - /i)^ + ^(fc)'*A2 v^(fc + /i)2 + F(fc)4A2 



(73) 



(74) 



From ( [7lD we see that the physical gap, namely the minimum energy of the 
quasiparticles, is F(/Lt)2|A|. Creating a quasiparticle-quasihole pair requires at 
least twice this energy. If we take F{p) to be a step function as above, then 
the gaps are simply |Ai| and lAgj. 

In this derivation, the variational nature of the procedure is somewhat 
hidden. The point is that the 9 angles were chosen to guarantee the absence 
of linear terms in the effective Hamiltonian. Such terms would indicate the 
instability of the no- "particle" state, where particles are defined by the action 
of the operators (^9|). This state varies with the 6 angles. Thus in searching 
for a ^-dependent Hamiltonian free of linear terms we are in effect looking for 
the stationary point of the expectaton value of the energy, within a continuous 
manifold of candidate states. 

For the same choice of parameters as we used in the 2SC case (namely 
A = 800 MeV, G chosen so that the vacuum constituent quark mass is M = 
400 MeV, and ji = 400 MeV) solving the coupled gap equations ( [73| ) yields 
Ag = 80 MeV and Ai ~ —176 MeV. As promised, this is not far from the ratio 
Ai — — 2A8 required if the condensate were entirely in the color 3 channel. 
The color 6 condensate is nonzero but small, as is also the case in the CFL 
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phase at asymptotically high densitiesE2lO If /i is increased to 500 MeV, we 
find Ag = 109 MeV and Ai = -249 MeV. These estimates of A may be a bit of 
an overestimate, because it is likely that the instanton interaction contributes 
significantly to M in vacuum but, because it is a six-fermion interaction, it 
does not contribute significantly to A. For this reason, we are overestimating 
G by assuming that M in vacuum is due entirely to the one-gluon exchange 
interaction. (This uncertainty did not arise with two flavors because in that 
case, the instanton and one-gluon exchange four-fermion interactions both con- 
tribute to M and to A in the same way.) Reducing G by a factor of two results 
in As = 16 MeV and Ai = -34 MeV at /i = 500 MeV. The magnitude of the 
gap is sensitive to the strength of the interaction, even though it is insensitive 
to details of the form factor as we described above. This reminds us that we 
should use these models as a qualitative guide only. 

In three-flavor QCD, the instanton interaction is a six-fermion interaction. 
As we have discussed in Section 2, its effects are small but are nevertheless 
important because they are the leading source of U{1)a breaking. It is also 
interesting to note that the four-fermion interaction introduced by instanton- 
anti- instanton pairs contributes to both M and A. In the analysis of Ref. 9, this 
interaction is employed in an analogous fashion to the way we have employed 
single-gluon exchange, and quite similar gaps are obtained. 

We can now outline the argument that the CFL phase is indeed the favored 
pattern of condensation. First, the interaction favors antisymmetry in color. 
Second, pairing with zero angular momentum is favored because it allows con- 
densation which utilizes the entire Fermi surface. This has been confirmed in a 
variety of explicit calculations. Pauli then requires flavor antisymmetry. These 
arguments favor a condensate of the form 

(^r^75V'f ) - e"^^6„-B ^ c^i . (75) 

Color-flavor locking corresponds to 0^ ~ diag(l, 1, 1), in which all nine quarks 
pair. The antithetical possibihty is cf)^ ~ diag(0,G, 1), but this is just the 
2SC phase, whose gap equation we have also solved! For ^ = 400 MeV and 
G fixed so that M = 400 MeV in vacuum, we found that in the CFL phase 
(with 4>'^ ~ diag(l, 1, 1)) there are 8 quasiparticles with gap 80 MeV and one 
with gap 176 MeV. For 0^ ~ diag(0,0, 1), on the other hand, we find four 
quasiparticles with gap 106 MeV and 5 quarks (strange quarks; blue up and 
blue down quarks) with zero gap. Clearly, the condensation energy (i.e. (fi) 
relative to that in the unpaired state) will be greater with 0^ ~ diag(l, 1, 1). 
Explicit evaluation of (^'|ri|\E') using the variational wave functions for the 2SC 
and CFL phases confirms this. 

These calculations are strongly suggestive that (f)^ ~ diag(l,l,l) is the 
choice with the lowest free energy, although they do also demonstrate that a 
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small admixture of the color 6 condensate which is symmetric in both color 
and flavor is mandatory. This means that a complete demonstration is more 
complicated. Schafer has, however, constructed an effective potential for </)^ 
(and the color 6 condensate) within a particular NJL model but making the 
most general assumptions for the form of the condensate. He finds that the 
CFL phase, with c/*;^ ~ diag(l, f, 1), minimizes (il). The authors of Refs. 62,12 
have addressed this question in analyses done using the interaction generated 
by the exchange of a propagating gluon, rather than in a model with point-like 
interactions. They confirm that (j)-^ ^ diag(l,l,l) is favored, thus explicitly 
confirming that the CFL phase is favored at arbitrarily high density. The 
essence of all of these different analyses is that in the CFL phase all quarks 
are able to pair, and the condensation energy is therefore larger than that for 
any less symmetric condensate. This simple argument could in principle be 
outweighed if a less symmetric condensate were found whose gaps were much 
larger in magnitude than those of the CFL phase. With all parameters fixed, 
though, different choices for (j)^ turn out to yield comparable gaps. Although 
(f)g ^ diag(0, 0, 1) does yield somewhat larger gaps, this advantage is not close 
to outweighing the 9/4 numerical advantage enjoyed by the CFL phase. 

4-3 Model Hamiltonian by Diagrammatic Methods 

In this brief subsection, we use diagrammatic methods, sometimes named after 
Dyson and Schwinger or, in this context, Nambu and Gorkov, to rederive the 
gap equation ( p2| ) for the 2SC phase, still in a model in which the QCD interac- 
tion between quarks has been replaced by a point-like four-fermion interaction. 
We do so both because these methods may be more familiar to some readers 
and because they serve as a warmup for the weak-coupling calculations of the 
next subsection. To that end, we also discuss the G ^ limit of the model. 

We follow the standard Nambu-Gorkov formalism and introduce an eight- 
component field {tp,'tp'^). In this basis, the inverse quark propagator takes the 
form 

where A = 7oA^7o and where in this expression A is a matrix with color, 
flavor and Dirac indices which have all been suppressed. The diagonal blocks 
correspond to ordinary propagation of massless quarks and the off-diagonal 
blocks reflect the possibility for "anomalous propagation" in the presence of a 
diquark condensate. The 2SC ansatz for the form of the gap matrix is 

Atf{p) = e'^''h.,Cj,A, (77) 

where A on the right-hand side is now the gap parameter, and all matrix 
structure has been written explicitly. 
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The full propagator of the interacting fermion system determines the proper 
self energy E via 

S^{S^'+Y.y^ (78) 

where S^^ is the inverse propagator in the absence of any interactions (given 
by ( [76| ) with A = 0) . The gap A is then determined by solving a self-consistent 
Dyson-Schwinger equation for S. Assuming a four-fermion interaction, this is 
given schematically by 




where the loop denotes a momentum integration over the full propagator S and 
where external legs have been amputated. Note that this equation sums an 
infinite class of diagrams, because S occurs within the propagator on the right- 
hand side. It is nevertheless (for obvious reasons) called the one-loop Dyson- 
Schwinger equation. As we shall see, it is precisely equivalent to the variational 
approximation we described in Section 4.2, in which we wrote an ansatz which 
admitted pairing but neglected other correlations. It is also often called a 
mean-field approximation. The Bogoliubov-Valatin derivation, also presented 
in Section 4.2 and also completely equivalent, makes this nomenclature clear. 

Using the four-fermion interaction ( |5^ ) of Section 4.2 which models single- 
gluon exchange, the Dyson-Schwinger equation is given explicitly by 

^ = -^G j -^aKS(j>)T^^ (79) 
where is the quark-gluon vertex 

^'^ " -(7mA^/2)^ ) ' ^^^^ 

and where we have chosen to work in Euclidean space. After some algebra 
(essentially the determination of S given the S^^ specified above), and upon 
suitable projection, this matrix equation reduces to a gap equation for the gap 
parameter A given by 

^^^^ j J2^\pl + {\p\-^lY + A^^ pl + {\p\ + ^i)^ + ^^) ' ^^^^ 

Whether the po integration contour is closed in the upper-half po-plane or in 
the lower-half plane, one pole from each of the two terms on the right-hand 
side contributes and one obtains the 2SC gap equation (^2|). The second term 
on the right-hand side of (Is^) describes the contribution of antiparticles to the 
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gap equation. The first term, due to particles and holes, yields a logarithmic 
divergence in the |p| integral if A ^ 0. As we have seen, this divergence at the 
Fermi surface is the hallmark of the BCS phenomenon. 

In the variational approach, once we have solved the gap equation we 
know the variational wave function and with this in hand we can evaluate the 
expectation value of the free energy If, instead, we have obtained the gap 
equation diagrammatically, we can still obtain the free energy upon realizing 
that the gap equation can be seen as the statement that the free energy is 
stationary with respect to the gap parameter. That is, we integrate the gap 
equation: 

/•Aaoiution / 2A f d^v \ 
n = rffrec + V J dAi^-— + 8j integrand j , (82) 

where the "integrand" is that on the right hand side of the gap equation (^l|), 
where Agoiution is the value of A which solves the gap equation, and where 
f2freo is that for noninteracting fermions. 

Our next goal is the study of color superconductivity in the limit of asymp- 
totically high density and thus, because of asymptotic freedom, at weak cou- 
pling. It goes against the spirit of the phenomenologically motivated models 
which we have used in Sections 4.2 and 4.3 to apply them at these densities. 
Our approach has been to use such models to investigate how the quarks re- 
solve the instability revealed by the renormalization group analysis of Section 
4.1 (they pair) and then to use the models as devices which allow us to make 
a phenomenological normalization of the strength of the attraction between 
quarks, and thus to obtain a qualitative measure of the magnitude of the gap. 
At asymptotically high densities, instead of attempting to normalize the in- 
teraction via vacuum phenomenology, we should use the weak-coupling QCD 
interaction itself. This is the subject of Section 4.4. Before beginning, and in 
order to provide a contrast, let us see what does happen if we take the four- 
fermion coupling G ~^ 0, and thus take A ^ /^,A. Solving the gap equation 
( |6^ ) explicitly in this limit yields 

A = 2V^»p(^),xp(-^) . (83) 

As expected from BCS theory, the factor 2/i^/7r^ which multiplies G in the 
exponent is the density of states (for red and green u and d quarks of a single 
chirality) at the Fermi surface in the absence of interactions. The A-dependent 
prefactor is clearly model-dependent. However, since the four-fermion inter- 
action with coupling G models single gluon exchange, the result ( ^3| ) could 
suggest that for /i ^ oo in QCD, we may find A ~ exp(— const /g^). We shall 
see in the next section that this is not correct. 



50 



Handbook of QCD / Volume 3 



The free energy can be obtained explicitly in the G ^ limit from ( p2| ) 
and, after using (p3[) to eliminate G, turns out to be 



V ~ 



It is nice that G and A, which pertain to a specific model, have disappeared. 
Only the physical observable A remains. It is reasonable to expect that in 
weak-coupling QCD, once A has been determined, the result ( ^4| ) is — except 
perhaps for the prefactor — a good estimate of the asymptotic condensation 
energy. Indeed, this can be argued for by dimensional analysis. In a similar 
spirit, we expect that the conclusions we have drawn in the previous section 
for what condensation pattern is favored in two- and three-flavor QCD are 
correct at asymptotic densities. Both these expectations are borne outt3 What 
remains to be determined is the asymptotic, weak-coupling, value of A. Given 
its model dependence, it is not surprising that the G ^ result ( |8^ ) for A 
fails to desctibe QCD at asymptotic densities; what was unexpected until the 
work of SoiO was that it fails qualitatively. 



4-4 Asymptotic Analysis 

Using diagrammatic methods, we now address the weak-coupling behavior of 
the gap A in QCD, which we expect will correctly reproduce the behavior of 
QCD at asymptotically high densities. We replace the point-like caricature 
of single-gluon exchange by single-gluon exchange itself. When retardation 
or relativistic effects are important, a Hamiltonian treatment (as in Section 
4.2) is no longer appropriate. One must pass to Lagrangian and diagrammatic 
methods, as in Section 4.3. The Dyson- Schwinger equation becomes 



A proper discussion of the fully microscopic calculatioiSiiffliiiSili 
is necessarily quite technical. Several subtleties arise, fundamentally because 
(as we have already seen in previous sections) ordinary perturbation theory is 
based on an unstable and qualitatively incorrect ground state and also because 
even if the ground state were valid there would be many charged low-energy 
excitations whose response produces electric screening and magnetic Landau 
damping, which drastically affect the nature of the interaction, and must be 
incorporated from the start. Before launching in, we present a schematic cal- 
culation, that conveys the spirit of the thing — and one of the most striking 
results — in easily digestible form. 
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Upon replacing gluon exchange by a contact interaction, we have already 
seen that we obtain a gap equation of the form 



where e = (|g| — /i) is the distance to the Fermi surface. The integral diverges 
at small e, so that as long as the proportionality constant is positive one will 
have non-trivial solutions for A, no matter how small is g. Indeed, for small g 
the solution is A ~ g"^"""*/"' . 

At asymptotically large /Lt, models with short-range interactions are bound 
to fail because the dominant interaction is that due-io the long-range mag- 
netic interaction coming from single gluon exchangewO When we now restore 
the gluon propagator, we find a nontrivial angular integral which diverges log- 
arithmically for forward scattering in the absence of any mechanism which 
screens the long-range magnetic interaction. We suspect, however, that in the 
superconducting phase this divergence will at the very least be precluded by 
the generation of a gluon mass oc A by the Meissner-Higgs mechanism.^ We 
therefore expect a gap equation of the form 



where is the angle between the external momentum and the loop momentum 
and where the integral is dominated by small e and small 6. This gap equation 
yields A - ^^^(log A)^ and hence A - e-'=°""'/9! 

The proper discussion of the microscopic gap equation is considerably more 
involved than (^6|). For example, it turns out that the coUinear logarithmic 
divergence is cutoff by dynamical screening (Landau damping) and not by the 
Meissner effectLJ As a result, the prefactor of 1/ff in the exponent turns out 
to differ from that which one would derive from ( pq ) . It nevertheless turns out 
that the conclusion that at weak coupling the gap goes, exponentially in the 
inverse coupling (rather than its square) still emergescS As a result, the gap 
is parametrically larger at — > oo than it would be for any point-like four- 
fermion interaction. This has the amusing consequence, that at asymptotically 
high densities the gap becomes arbitrarily large! This is because asymptotic 
freedom ensures that it is the microscopic coupling l/g(/i)^ which vanishes 
logarithmically, so that e"'^"'^''*/^^'^) does not shrink as fast as Since the 

"dimensional analysis" scale of the gap is set by /i, A grows asymptotically 
even though A//z shrinks. 

The Meissner mass is of order g/i but the gluon wave function renormalization 

is of order g/^/'A and the physical screening length, which is the ratio of the two, is therefore 
of order l/AEj 




(85) 




(86) 
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We shall discuss the asymptotic analysis by example, focussing on the 
2SC phase and following the presentation of Ref. 64. After describing the 
results, we shall quote the analogous results for the CFL phase. We use the 
2SC phase as our example only because the algebra is simpler with only two 
flavors. The physics is under much better control in the CFL phase. In the 
calculation we outline, we determine the dominant gap in the 2SC phase, 
describing the pairing of red and green u and d quarks. This leaves the blue and 
strange quarks unpaired, and thus leaves open the possibility of further infrared 
divergences at the Fermi surface, leading to small gaps via the formation of 
Cooper pairs with nonzero angular momentum. Also, three gluons are left 
unscreened raising the spectre of further infrared complications. In the CFL 
phase, the dominant condensate gives a gap to all nine quarks and gives a mass 
to all eight gluons, and there are no further infrared difficulties to be resolved. 



As in Section 4.3, we introduce an eight-component field (■0, w") and 
the inverse quark propagator ([76|). We assume that the gap matrix is anti- 
symmetric in both color and flavor, which is the channel in which single-gluon 
exchange is attractive. We also assume that the condensation is in the channel 
with total angular momentum J = 0. In the case of short range interactions, 
all these assumptions can be justified by reference, to the renormalization group 
equations for a general four-fermion interactionpHH as described in Section 4.1. 
The single-gluon exchange interaction is long range, however, and other forms 
of pairing might take place in addition. In particular, since the interaction 
is dominated by almost coUinear scattering, it includes substantial strength 
in all partial pSsiajies, and one might expect some condensation in the higher 
partial wavesEBEil It turns out that theae-cpiadensates have the same leading 



g-dependence as the J = condensatepj'LitEj but they are suppressed by nu 
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merical factors of order 1/100 to l/lOOoM We therefore neglect them, and 
concentrate on the J = gap. 



We also assume that the gap has positive parity. As we discussed in 
Sections 4.1 and 4.2, single-gluon exchange does not distinguish between the 
positive parity (tpC^stp) condensate and the negative parity (tpCtp) conden- 
sate. This Jfigeneracy is lifted by instantons, which favor the positive parity 
channelui3'E3. At large chemical potential instanton effects are exponentially 
suppressed. In the following, we will therefore assume that the only instanton 
effect is to determine the parity of the gap. Finally, we neglect quark mass 
effects and chiral symmetry breaking LR condensates. As shown in Ref. 58 
and described in Section 4.1, there is no BCS instability in the case of pairing 
between left and right handed, quarks. The formation of LR condensates is 
therefore suppressed by m//i Ell 
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.all these assumptions and specializations, the gap matrix takes the 



A,f b) = (A2)°''(t2).,C75 (Ai(po)P+(p) + A2bo)P-b)) , (87) 
where the projection operators are given by 

P+(p) = i(l+7o7-p), 

P.{p) = i(l-7o7-p). (88) 

We have neglected the dependence of the gap on the magnitude of the mo- 
mentum IpI, but kept the dependence on frequency po- The |p|-dependence 
can be dropped because, in the weak coupling limit, all momenta are close to 
the Fermi surface. For short-range interactions, the dependence on frequency 
can also be neglected. Thus, in Section 4.3 we chose the ansatz ( |77| ) in which 
Ai(po) — ^2[po) and both were given by the constant A. This was appropriate 
in that context, and we had no difficulty finding a solution to the gap equation 
upon making this ansatz. Now, however, because long range interactions are 
important, retardation effects cannot be neglected. 

P+ and P- are projection operators which project onto particles and an- 
tiparticles respectively. This means that Ai(po) describes the modification of 
the propagator due to particle-particle and hole-hole pairing whereas A2(po) 
describes that due to antiparticle-antiparticle pairing. The gap in the quasi- 
particle spectrum at the Fermi surface is given by Ai(Ai), and we are therefore 
primarily interested in deriving a gap equation for Ai(po)- We expect that at 
weak coupling, quarks near the Fermi surface should dominate and the contri- 
bution of A2 to the gap equation for Ai should vanish. (This is what we found 
in our analysis of the gap equation (|6^), for example.) This emerges explicitly 
as a consequence of the fact that in the weak coupling limit, we can replace 
7o7 • p by the unit matrix using the equations of motion. Note that this does 
not mean that A2 is in any sense smaller than Ai. It only means that the 
contribution of A2 to the gap equation for Ai becomes small. 

The self energy in the Nambu-Gorkov formalism obeys the Dyson-Schwinger 
equationQ 

= J -^Tp{p)TtD,,{k - p) . (89) 

Here, S(fc) is the proper self energy defined via ([78|), is the quark-gluon 
vertex and D^,y{k — p) is the gluon propagator. We have written ( |89| ) in 
Minkowski space. To leading order in the perturbative expansion, we can use 
the free vertex (^0|). To leading order, we can also neglect the diagonal part 
of the proper self energy, namely the fermion wave function renormalizationEj 
We shall see below that although we obtain the asymptotic form for the gap, 
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the results we obtain are uncontrpllcd for g > gc ^ 0.8, corresponding to a 
rather large value of /x ~ 10* MeVEJ This breakdown could in principle reflect 
a failure in any of our assumptions. We expect, however, that it arises because 
contributions which have been truncated in writing the one-loop Schwinger- 
Dyson equation ( |89| ) with (|80| ) are large for g > gc- That is, we expect that 
this truncation (for example the neglect of vertex corrections) , and not any of 
the simplifications we introduced in making the ansatz (|8^) for A, is the most 
significant assumption we make in this analysis. 

The Dyson-Schwinger equation ( p9| ) reduces to an equation for the gap 
matrix. 



S2i{p){i,— ]D^,{k-p). (90) 



Here, 5*21 (p) is the 21-component of the fermion propagator in the Nambu- 
Gorkov representation. 5*21(^1) is determined from the inverse of (^. We 
have 

1 1_ 

Inserting the ansatz (^) for the gap gives 



S2l{p)^- ,^ .,.T ^7T-^. Xu^, .nri-iA - (91) 



S2i{p) = -A2-r2C75 



Ai(po)P-(p) , A2{po)P+(p) 



.pI " m - y^f - Ai(po)2 pI - m + - A2(po)2 

(92) 

Both the RHS and the LHS of the gap equation are proportional to T2, so the 
flavor structure simply drops out. The color coefficient is given by 

^(Aa)^A2A, = -^^A2 = -^A2 (iV, = 3), (93) 

where we have used the Fierz identity (A")ij (A")*;; = —{2/Nc)SijSki + ^SgSjk 
and the factor 1/4 comes from the color generators = A°/2. Projecting (90) 
on Ai_2 gives two coupled gap equations 



hr (7^P+(p)7,P^(fc)) 2 n^^^t'^ a ( ^2 \ D^,{k-p) , (94) 
2 Po - \\P\ + N - A2(po)^ 



where the upper and lower signs on the RHS correspond to Ai and A2 on the 
LHS. 
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We now must specify the gluon propagator. The giuon propagator in a 
general covariant gauge is given by 



DM = + - , (95) 

q - G{q) r - F[q) 9* 

where the projectors Pj^^^ are defined by 

PT^=5,,-m,, Po^o=Po^ = 0, P^^^ ^ -g^, + ^ - Pj^ . (96) 

The functions F{q) and G{q) describe the effects of the medium on giuon prop- 
agation. If we neglect the Meissner effect (that is if we neglect the modification 
of F{q) and G{q) due to the gap A in the fermion propagator) then F{q) de- 
scribes Thomas-Fermi screening, G{q) describes dynamical screening (Lpadau 
damping), and they are given in the hard dense loop approximation by til 

Fiq) . 2.n^^|l-^QofelV Qo(.) = ^ log fi^i^ V (97) 




G{q) = m — 



\q\^\\q\ r 2"^\x-i 



i HiM 



(98) 



where rn^ — N f g'^ fi^ / {An'^) and we are working with Nf — 2 flavors. Note 
that G ^ for (7o ^ 0: static magnetic modes are not screened. Nonstatic 
modes are dynamically screened due to Landau damping. Note that the gluon 
propagator contains the gauge parameter ^, which must not appear in physical 
results. 

In the weak coupling limit, go is small as compared to \q\. In this case we 
can expand the projectors P^^ ~ Sf^oS^o and qiqj/q'^ — qiqj- (See Ref. 69 for 
an analysis in which these simplifying kinematic approximations and others 
which follow are not made.) The gap equation now becomes 



Ai(fco) 



2z.g2 r d^p j Ai(po) 



i27r)^\pl-m-f^r~^i{poy 



I ■ p \ + \k-p 



{k-pf-G{k-p) [k-pY -F{k-p) 



A2(Po) 



Pl-m+^^?-^2{J>,y 



\ + jk-p J - ^k -p ^ 



[k - pf - G{k - p) [k - pY - F{k - p) [k - pf 



(99) 
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There is a similar equation for A2 in which the two terms in the round brackets 
are interchanged. Only the first term in ( p9| ) has a singularity on the Fermi 
surface. In the weak coupling limit, we can therefore drop the second term, 
and we are left with an equation for A(A:o) = Ai(fco)- This equation is in- 
dependent of the gauge parameter ^. The ^-independence of the Ai term in 
( p9| ) arises because of the kinematic approximations made_in deriving it. If 
one goes beyond these approximations, Ai is ^-dependents The second gap 
parameter A2 is not suppressed in magnitude and its gauge-dependence is not 
kinematically suppressed. However, A2 does not lead to a gap on the Fermi 
surface. 

The fact that the gap is gauge independent in the present weak-coupling 
approximation is a consequence of the fact that the gap is determined by the 
scattering of quarks that are almost on shell. For on-shell quarks, the fact that 
the gauge dependent part of the propagator does not contribute follows directly 
from the equations of motion for the quark fields. This argument is adequate 
only to lowest order, as appeared in the specific diagrammatic cancellation 
(and non-cancellation) above. 

For large chemical potential the integral over p is dominated by momenta 
in the vicinity of the Fermi surface, |p| ~ /i and pq <^ 11. We can expand all 
momenta as p = pp + I, where pp is on the Fermi surface, and I is orthog- 
onal to it. Asymptotically, |/| <C \pf\ and the integration measure becomes 
dpo ji^ did cos 9 d(j). We also have |A: — p| ~ •\/2/x(l — cosO). The integral over 
is performed trivially. We analytically continue to imaginary p^, and perform 
the integral over / by picking up the pole in the diquark propagator. We find 



A(fco) = j dpa J d cos 61 



I - i cos 6* 



l-cos9+{G+{ko-poy)/{2^i^) 
i + ^cos0 ^ A{po 



1 - cos H- (F + {ko - Po)2)/(2Ai2) J + A(po) 



(100) 



The integral over cosO is dominated by small 0, corresponding to almost 
collinear scattering. We must, therefore, take careful account of the medium 
modifications of the gluon propagator described by F{q) and G(g) given in (p7|). 
In the gap equation ( |lOO| ) F{q) and G{q) are to be evaluated at go = — Vo 
and \q\ = \k — p\ ~ \/2/^(l — cos 9), and are therefore functions of fco — pq and 
cos9. For go ^ 9* ^ the expressions ( ^t]) for F{q) and G{q) simplify, yielding 

F{q) = 2m\ G{q)^^m^^. (101) 

In the longitudinal part, m\, = 2m? is the familiar Debye screening mass. In 
the transverse part, nonstatic modes are dynamically screened. In our case. 
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typical frequencies are on the order of the gap, go — This means that the 

1/2 

electric part of the interaction is screened at qs — rrij^ whereas the magnetic 
interaction is screened at qm — {7r/4 • m|jA)^/'^. 

Asymptotically, qm Qe, and magnetic gluon exchange dominates over 
electric gluon exchange. We therefore begin by analyzing the gap equation 
taking into account the magnetic part of the interaction only. We will also 
approximate cos9 ~ 1 in the denominator and drop (/cq — Po)'^ in the denomi- 
nator. All of these terms will be reinstated later. The integration over cos 6 is 
now straightforward. We have 

IStt^ J \ Nfg^\ko-po\J y/pl + A{pn)^ 

If we are only interested in the leading exponential behavior of the gap we 
can drop the numerical factors and the powers of g in the logarithm. We then 
arrive at 

ISlT^ J V 1^0 -Pol/ y/pl + A{po)^ 

where c is a factar to be determined more fully below. This equation was first 
derived by Sonu2l (In ordinary|_s.uperconductivity, the corresponding equation 
was first derived by EliashberglS) What we have shown here (following Ref. 64) 
is that one can derive this equation directly from the Dyson-Schwinger equation 
in the weak coupling limit, and that upon making the approximations we have 
made the result is independent of the gauge parameter. The integral equation 



(103) can be converted to a differential mualiipnEj and in the weak coupling 




A(fco) ~ Ao sin (-^ log (^) ) , fco > Ao , (104) 



limit an approximate solution is given b; 

where ^ 

Ao = 2cfiexp (^"-^^ • (105) 

At this point a few comments are in order. First, we note that the use 
of perturbation theory to determine the dynamic screening is self-consistent. 
Since A ~ /iexp(— constant/g), the gap grows as /z ^ c» and qM ^ ^qcd- 
Second, we note that it is essential to keep the frequency dependence of the gap. 
For small frequencies, A{ko) varies over scales on the order of kg ~ Ao itself. 
Therefore, A(A;o) cannot be replaced by a constant. Were we to approximate 
A(fco) — Ao, we would obtain a gap equation for Aq that has the correct double 
logarithmic structure and gives Ao ~ /iexp(— constant/t/), but the constant in 
the exponent would not be correct. 
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We now come to the role of electric gluon exchanges. We include the 
second term in (IOC) with F = m|,. We again use the approximation cos 9 ~ 1 



in the numerator and drop the (fco~Po)^ term in the denominator. Let us note 
that in the forward direction, electric and magnetic gluon exchanges have the 
same overall factor. Performing the integral over cosO, we find 

^ ^/dp„{iogfi+ „,ir,., )aiogfi 



(106) 



where the factor 3/2 in front of the second term comes from the difference 
between dynamic screening, ~ l?!^''^, an d st atic screening, qs ^ \q\- In 



the weak coupling limit, the gap equation (106) reduces to the form (103) 



with approximate solution (104), where we now see that we can estimate the 
constant 2c to be 

2c = 2048%/27r'*iV^^^/^ = hlin'^g-^ ~ 5.0 x 10'^ g'^ {Nf = 2) . (107) 



We have obtained the result (107) for the prefactor 2c in the expression for 
the gap Aq by collecting the leading logarithms from both electric and magnetic 
gluon exchange. To do better, Schafer and Wilczek solve the gap equatiati 
( |100| ), which incorporates both electric and magnetic exchange, numerically^ 
They also keep the cosO dependence in the numerator, and the terms (fep— po)^ 
in the denominator. Finally, they use the exact forms of G{q) and F[q) in the 
hard dense loop approximation, (^7|). This takes into account that there is 
no dynamic screening for |g| < q^. Asymptotically (i.e. for small g) the 
numerical solution to the gap equation is well described by (|104|105D with 



2c ~ 1.4 X 10^5-5, at least for ko < y/Nj /{S'K)gfi. At lar ger ko, the retardation 
terms ~ (fcp — Po)^ dominate over screening and A falls off more quickly with 



increasing fcg than in (104). We are interested in ko — Aq, however, and in 
this regime (104) describes the shape of the numerical solution quite well. 

In order to use the weak-coupling results, we must choose an energy scale 
at which to evaluate the running QCD coupling. Schafer and Wilczek chose to 
use the one-loop running coupling constant evaluated at the Fermi momentum 
Pp — fi, which is an average over the momenta of the exchanged gluons, which 
are in the range [qM , 2/i] . Strictly speaking, only a higher order calculation 
can fix the scale in the running coupling. It turns out, however, that with g 
evaluated at /i, the gap Ag only depends rather weakly on /i. We return to 
this below. The numerical solution yields Aq ~ 40 MeV at /i = 10^" MeV, 
corresponding to Aq/m = 4 x 10~^ at g = 0.67 and a value of 2c which is 



about 2/5 that in (107), and Ao ~ 90 MeV at /i = 400 MeV, corresponding to 
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Ao//i = 0.23 at 5 = 3.43 and a value of 2c which happens to be almost exactly 
that in ( |107D . As ^ increases from 400 MeV to about lO*" MeV, A drops from 
90 MeV to about 10 MeV; as ^ increases further, A then begins its inexorable, 
but logarithmic, rise. 

Clearly, with /i = 400 MeV and g = 3.43, the calculation of the gap is 
determined by momenta which are not large compared to Aqcd, and all of 
the approximations we have used have broken down. More on this below, but 
before plunging into caveats let us note that it is gratifying to see that the 
order of magnitude of the result agrees with that obtained in the calculations 
described in Sections 4.2 and 4.3, based on more phenomenological effective 
interactions, which were normalized to the strength of chiral symmetry break- 
ing at zero density, rather than to the calculable asymptotics of the running 
coupling. _ 

Son's result e3 




(108) 



has now been confirmed using a variety of methods 0Elll'0ilil'0iiii How- 
ever, even the C((7°) contribution to the prefactor b in ( |l08| ) is not yet fully un- 
derstood, and certainly nobody has proposed a controlled calculational scheme 
which could be used to push the calculation of b to arbitrarily high accuracy. 
We have worked through the analysis of Schafer and WilczekEJ which leads to 
the estimate that b ^ 5127r^ in the 2SC phase, a result yd|U£|l|^pS|Upp by 
their numerical work and by the work of other authorsEil'Ea'E3ll3C3 A similar 
analysis yields the estimate that b ~ 5127r''2~^/^(2/3)^/^ in the CFL phase, 
with the (2/3)5/2 

coming from the iV/-dependence of the Debye mass and the 
2~^/^ copaing from an analysis of the color- flavor structure of the CFL gap 
equationt3 along the lines of that we have presented in Section 4.2 in the 
context of models with four-fermion interactions. Numerical solutions of the 
weak-coupling CFL gap equation by Evans et al. are in good agreement with 
this estimate, and these authors also confirm that including Meissner effects 
in the gluon propagator makes little difference to the final resultE3 

The argument that the CFL phase is favored at asymptotically high density 
in QCD with three flavors of quarks is similar to that given at the end of Section 
4.2. 2SC pairing would yield a gap given by (|lO|) with b ~ 5127r4(2/3)5/2, 
whereas in the CFL phase the gap is smaller by a factor of 2^^^^. In the 
asymptotic regime, the_miidensation energy is proportional to A^/i^ up to 
logarithmic correctionsllj'tj As expected, the expression ( |8^ ) derived by taking 
G — > in the analysis done using a four-fermion interaction is a good guide. 
In the 2SC phase, four quarks pair with gap A and the condensation energy 
energy is ~ 4:{A'^ ^'^ / Att'^) . In the CFL phase at asymptotically high density, 
the color 6 condensateJs suppressed relative to the color 3 condensate by a 
factor f/'\/21og(2)/367rliHl and so to leading order one has eight quarks with 
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gap A and one with gap 2 A, and a condensation energy which is therefore 
~ 12(A2/x2/47r2). Although the CFL gap is smaUer than the 2SC gap by a 
factor of 2~^/'^, this does not outweigh the advantage the CFL phase enjoys 
because it ahows all nine quarks to pair. As a result, the CFL phase is favored. 
As we saw in Section 3, this conclusion continues to hold even when the strange 
quark mass is nonzero, as long as m^/2/z is small compared to A. 

There are several contributions to b which we have not taken into account 
above. They modify the magnitude of the gap, but do so in the same way for 
CFL or 2SC pairing. For example, modifications to the quasiparticle dispexsipn 
relations in the normal (nonsuperconducting; high temperature) nhaseEaSa 
and quasiparticle damping effects in the superconducting phase 113 do affect 
the gap, and both tend to reduce b. In the normal phase, the effects of wave 
functiop-renormalization are only important within ~ exp(— l/t;^) of the Fermi 
surfacel^ which supports Son's conclusion that in the^aiiperconducting phase 
these effects result only in a subleading CQjtrection to AEj This correction does, 
however, reduce b by about a factor of 5E3 

The value of b is also affected by the choice of the scale at which g is 



evaluated in (W8). The consequent modifications to A have been considered 
by Beane et a/H They do a rjeaormalization group analysis within the effective 
theory constructed by Hongp3 who noted that the antiparticles can safely be 
integrated out. Beane et al use the renormalization group to compute and 
resum contributions of the form q;"+^/?q log"(A//z) to the gap equation, where 
/3o is the first coefficient in the QCD /3-function. Their results demonstrate 
that g should be evaluated at a /z-dependent scale which is much lower than 
/iEa Upon assuming that the QCD coupling runs according to the vacuum /3- 
function all the way down to the magnetic scale qM — an assumption which 
is not justified but which has not yet been superseded — they find that if, by 
convention, g is taken as ^(/x), then b is enhanced by a factor of about 20. 

Finally, we return to the question of the gauge dependence of b. Both 
the analysis of Schafer and Wilczek presented above and that of Pisarski and 
RischketJ demonstrate that the 0{g^) contribution to b is gauge- independent. 
However, Rajagopal and Shuster have gone back to (|9j), eliminated A2 in 
the gap equation for Ai, but made no further kinematic simplifications. The 
resulting gap equation for Ai is gauge-dependent. Examination of the gauge- 
dependent (and (7-dependent) contributions to b obtained by solving this equa- 
tion reveals that these do decrease with decreasing g, but that they only begin 
to decrease for g < 0.8E3 This means that effects like vertex corrections which 
have to date been neglected in all calculations based on the one-loop Schwinger- 
Dyson equation (e.g. those of Refs. 61,63,64,66) are small corrections to b only 
for n » 10^ MeV. 

All these caveats should by now have made it clear that although we are 
confident of the leading g-dependence of the gap A at asymptotically large 
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/i, there remains physics which has not been satisfactorily treated which con- 
tributes at order to the prcfactor h. Furthermore, effects which are higher 
order in g are large enough that they can only be neglected for /i ^ 10* MeV. 

The asymptotic results nevertheless have very important consequences. 
We have seen in Section 3 that the outcome of the competition between the 
CFL and 2SC phases depends on the relativ e ma gnitudes of A and m1/2iJL. The 
latter decreases at large /x, while the result (108) demonstrates that A increases 
logarithmically as /i oo. This means ihat the CFL phase is favored over the 
2SC phase for /i ^ oo for any TOs ^ oot3 Even though the asymptotic regime 
where A can be calculated from first principles with confidence is not accessed 
in nature, it is of great theoretical interest. The weak-coupling calculation of 
the gap in the CFL phase is the first step toward the weak-coupling calculation 
of other properties of this phase, in which chiral symmetry is broken and the 
spectrum of excitations is as in a confined phase. As we have described in 
Section 2, for example, the masses and decay constants of the pseudoscalar 
mesons can be calculated from first principles once A is known. 

Although the value of A is under control asymptotically, it seems fair to say 
that applying these asymptotic results at /j, = 400 MeV is currently at least as 
uncertain a proposition as applying estimates made using phenomenologically 
normalized models with point-like interactions. Nevertheless, if we take the 
estimates for the prefactor b provided by Schafer and Wilczek's numerical 
results described above and apply them at /i ~ 400 MeV, they predict gaps of 
order 100 MeV. The consequent critical temperatures are related to the ^tQ 
temperature gap A by the standard weak-coupling BCS result Tj. = O.STAjElEj 
and are therefore of order 50 MeV. We have seen that some known corrections 
push this estimate up while others push it down, and that the calculation 
whence it came is, regardless, of quantitative value only for ^ ^ 10* MeV. 
It is nevertheless satisfying that two very different approaches, one using zero 
density phenomenology to normalize models, the other using weak-coupling 
methods valid at asymptotically high density, yield predictions for the gaps and 
critical temperatures at accessible densities which are in qualitative agreement. 
Neither can be trusted quantitatively for quark number chemical potentials 
fjL ~ 400 — 500 MeV, as appropriate for the quark matter which may occur in 
compact stars. Still, both methods agree that the gaps at the Fermi surface 
are of order tens to 100 MeV, with critical temperatures about half as large. 

Tc ~ 50 MeV is much larger relative to the Fermi momentum (say /i ^ 
400 — 500 MeV) than in low temperature superconductivity in metals. This 
reflects the fact that color superconductivity is induced by an attraction due 
to the primary, strong, interaction in the theory, rather than having to rely 
on much weaker secondary interactions, as in phonon mediated superconduc- 
tivity in metals. Quark matter is a high- Tc superconductor by any reasonable 
definition. It is unfortunate that its Tc is nevertheless low enough that it is 
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unlikely the phenomenon can be realized in heavy ion collisions. 
4-5 Challenges for the Future 

While the rigorous asymptotic justification of color superconductivity is a mar- 
velous result, there are compelling physical and mathematical motivations to 
go further. 

For physical applications, wc arc interested in densities not much beyond 
nuclear, which (as we have seen) appears to be quite far from asymptopia. It 
is not at all clear, in particular, that the form of the interaction which dom- 
inates asymptotically, namely near-forward magnetic scattering, is the most 
important at physically achievable densities. Indeed, electric scattering — or 
binding through electric flux tubes! — and instantons more nearly resemble 
the interactions used successfully in phenomenological models of hadrons at 
zero density and applied to nonzero densities as described in Sections 4.2 and 
4.3. 

On the mathematical side, it clearly would be very desirable to have a sys- 
tematic approximation framework that did not depend on selective diagram 
resummation, and whose results were manifestly gauge invariant. Similar prob- 
lems arise in the treatment of ordinary superconductors in condensed matter 
physics, albeit in a simpler form since the gauge group is abelian. However, 
in that context the lowest-order approximation is adequate in practice since 
elcctrodynamic radiative corrections are overwhelmed, quantitatively, by many 
other uncertainties. 

In formulating an appropriate perturbation theory, the difficulties arise due 
to the gauge dependence of the primary order parameter and to the compli- 
cated (in particular, non-instantaneous) nature of the microscopic interaction. 
The whole notion of the primary order parameter depends on fixing a gauge, 
and we assume that is done. On physical grounds, as we have emphasized 
already, one expects a weak-coupling approach to work only when it is based 
on perturbing around the correct, gapped ground state; but of course this begs 
the question of how to construct the groimd state. The BCS-inspired varia- 
tional procedure works well for instanteous interactions, where one can define 
a simple Hamiltonian, but is very awkward for retarded interactions. As we 
saw in Section 4.2, in one of its many implementations (that due originally to 
Bogoliubov and Valatin) it takes the form of a search among "trial Hamiltoni- 
ans" incorporating quadratic terms mixing particles and holes, with coefficients 
depending on the momentum. This suggests that one should employ trial La- 
grangians with similar terms, allowing also for energy dependence, using the 
action variational principle popularized by Feynman in connection with the 
polaron problem. Such a procedure has the great advantage that one need 
only do perturbation theory around the ground state of the trial Lagrangian, 
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which can be designed to incorporate the correct physics (i.e. gaps) to begin 
with. In principle, the choice of gauge can simply be treated as an additional 
variational parameter. 

Ultimately, truly quantitative results on QCD at moderately ultra-nuclear 
densities will probably have to await numerical work using lattice gauge theory, 
as has been the case for zero density. In principle, we can bypass all the dif- 
ficulties involved in defining the primary condensate by working directly with 
the distinctive gauge-invariant parameters predicted for color-flavor locking. 
Specifically, one should encounter diverging susceptibilities for chiral symme- 
try breaking, especially in the diquark-antidiquark channel (^, and for U{\)b 
breaking, in the baryon number 2 channel. The corresponding vacuum ex- 
pectation values, signaling chiral symmetry breaking and superfluidity, could 
be stabilized by adding the corresponding inflnitesimal source terms. These 
features, and the associated appearance of Nambu-Goldstone bosons, would 
serve as unmistakable signatures of the color-flavor locked phase. In practice, 
unfortunately, all presently known numerical algorithms for QCD at flnite den- 
sity and zero temperature converge too slowly to be useful. The discovery of 
suitable algorithms for this problem stands as perhaps the greatest theoretical 
challenge in QCD today. 

5 The Phase Diagram 

In this article we focus largely on zero temperature quark matter, but it is 
important to understand where the 2SC and CFL phases flt on the phase 
diagram of QCD as a function of both temperature and density. Because QCD 
is asymptotically free, its high temperature and high baryon density phases 
are more simply and more appropriately described in terms of quarks and 
gluons as degrees of freedom, rather than hadrons. At high temperatures, in 
the resulting quark-gluon plasma (QGP) phase all of the symmetries of the 
QCD Lagrangian are unbroken and the excitations have the quantum numbers 
of quarks and gluons. At high densities, on the other hand, we have seen 
that quarks form Cooper pairs and new condensates develop. At high enough 
density, chiral symmetry is broken by color-flavor locking. At densities which 
are high enough that nucleons overlap and the matter is in a quark matter 
phase but which are not high enough for color-flavor locking, we expect to flnd 
the 2SC phase. 

In this section, we describe the qualitative features of the QCD phase 
diagram. We choose to describe the entire phase diagram, and not just the 
region of low temperatures and high densities. The exploration of the higher 
temperature regions of the diagram is the object of extensive experimental 
efforts in heavy ion collision experiments at CERN and Brookhaven. Due 
to its topical interest, we describe the high temperature regime of the phase 
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Figure 1: QCD Phase diagram for two massless quarks. Chiral symmetry is broken in 
the hadronic phase and is restored elsewhere in the diagram. The chiral phase transition 
changes from second to first order at a tricritical point. The phase at high density and 
low temperature is a color superconductor in which up and down quarks with two out of 
three colors pair and form a condensate. The transition between this 2SC phase and the 
QGP phase is likely first order. The transition on the horizontal axis between the hadronic 
and 2SC phases is first order. The transition between a nuclear matter "liquid" and a gas 
of individual nucleons is also marked. At T = 0, it separates the vacuum phase from the 
nuclear matter phase; Lorentz-boost symmetry is broken to its right but unbroken to its 
left. At nonzero temperature, Lorentz-boost symmetry is broken in both the nuclear gas 
and nuclear liquid, and this line of phase transitions may therefore end. It is thought to end 
at a critical point at a temperature of order 10 MeV, characteristic of the forces which bind 
nucleons into nuclei. 

diagram in some detail, in addition to explaining how the dense quark matter 
phases which occupy us in the remainder of this article fit in. In Section 6 
below, we refocus on cold dense quark matter, as we describe current efforts 
to understand how to use observed phenomena occuring in compact stars to 
map this region of the phase diagram. 

Let us begin with a brief review of the phase changes wich occur as a 
function of temperature at zero baryon number densityEJ That is, we begin by 
restricting ourselves to the vertical axis in Figures 1 through 4. This slice of the 
phase diagram was explored throughout the early universe during the first tens 
of microseconds after the big bang. It can also be studied in lattice simulations. 
As heavy ion collisions are performed at higher and higher energies, they create 
plasmas with a lower and lower baryon number to entropy ratio and therefore 
explore regions of the phase diagram closer and closer to the vertical axis. 

In QCD with two massless quarks {mu,d = 0; iris = oo; Figure 1) the 
vacuum phase, with hadrons as excitations, is characterized by a chiral con- 
densate {'tpLfpR)- Whereas the QCD Lagrangian is invariant under separate 
global flavor rotations of the left-handed and right-handed quarks, the pres- 
ence of the chiral condensate spontaneously breaks SU{2)]^ x SU{2){i to the 
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subgroup S'C/(2)i_|_/f , in which only simultaneous flavor rotations of L and R 
quarks are allowed. Locking left- and right-handed rotations in this way breaks 
global symmetries and results in three massless Goldstone bosons, the pions. 
The chiral order parameter, a 2 x 2 matrix M"'' in flavor space, can be written 
in terms of four real flelds a and tt as 

{i'lc.rii) = = aS^' + 7? . {rf , (109) 

where the f are the three Pauli matrices. SU{2)l and SU{2)ii rotations act 
on M'-' from the left and right, respectively. The order parameter can also be 
written as a four component scalar field (j> — {a, v?) and the SU(2)]^ x SU{2)^ 
rotations are then simply 0(4) rotations of cf). In this language, the symmetry 
breaking pattern S'[/(2)l x SU{2)r SU{2)l+r is described as 0(4) 0(3): 
in the vacuum, (0) ^ and this condensate picks a direction in 0(4)-space. 
The direction in which the condensate points is conventionally taken to be the 
a direction. In the presence of (a) ^ 0, the if excitations are excitations of the 
direction in which {(p) is pointing, and are therefore massless goldstone modes. 

At nonzero but low temperature, one finds a gas of pions, the analogue 
of a gas of spin waves, but {(f)) is still nonzero. Above some temperature Tc, 
entropy wins over order (the direction in which 4> points is scrambled) and 
((/)) = 0. The phase transition at which chiral symmetry is restored is likely 
second ordet,and belongs to the universality class of 0(4) spin models in three 
dimensionsE3 Below Tc, chiral symmetry is broken and there are three massless 
pions. At T = Tc, there are four massless degrees of freedom: the pions and the 
sigma. Above T = Tc, the pion and sigma correlation lengths are degenerate 
and finite. 

In nature, the light quarks are not massless. Because of this explicit chiral 
symmetry breaking, the second order phase transition is replaced by an ana- 
lytical crossover: physics changes dramatically but smoothly in the crossover 
region, and no correlation length diverges. Thus, in Figure 2, there is no sharp 
boundary on the vertical axis separating the low temperature hadronic world 
from the high temperaturejmark-gluon plasma. This picture is consistent with 
present lattice simulationsPE^ which suggest 21^ MQ-p 190 MevNa 

Arguments based on a variety of modcliEj'LiiLrO'lij'LJ indicate that the chi- 
ral symmetry restoration transition is flrst order at large fi. The 2SC quark 
matter phase that arises in two-flavor QCD at values of fi above this first or- 
der transition features new condensates, different from those in the hadronic 
phase, and these condensates do not break chiral symmetry. The fact that this 
is a transition in which two different condensates compete iitmngthens previ- 
ous model-based arguments that this transition is first orderOcEl This suggests 
that the phase diagram features a critical point E at which the line of first or- 
der phase transitions present for /i > he ends, as shown in Figure 2.^ At he, 

'^If the up and down quarks were massless, E would be a tricritical pointgj at which the first 
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Figure 2: QCD phase diagram for two light quarks. Qualitatively as in Figure 1, except 
that the introduction of light quark masses turns the second order phase transition into a 
smooth crossover. The tricritical point becomes the critical endpoint E, which can be found 
in heavy ion collision experiments. 

the phase transition is second order and is in the Ising universahty class@El 
Although the pions remain massive, the correlation length in the a channel 
diverges due to universal long wavelength fluctuations of the order parameter. 
This results in characteristic signatures, analogues of critical opalescence in 
the sense that they are unique to coUisiapa which freeze out near the critical 
point, which can be used to discover i?EJe3 

Returning to the ^ = Q axis, universal arguments^ again backed by lattice 
simulationpZI tell us that if the strange quark were as light as the up and down 
quarks, the transition would be first order, rather than a smooth crossover. 
This means that if one could dial the strange quark mass TOs, one would find 
a critical ,m1 at which the transition as a function of temperature is second 
orderEaO Figures 2, 3 and 4 are drawn for a sequence of decreasing strange 
quark masses. Somewhere between Figures 3 and 4, nis is decreased below 
and the transition on the vertical axis becomes first order. The value of 
is an open question, but lattirjp .simulations suggest that it is about Jialf the 
physical strange quark massE3t3 These results are not yet conclusive^ but if 
they are correct then the phase diagram in nature is as shown in Figure 3, and 
the phase transition at low /i is a smooth crossover. 

These observations fit together in a simple and elegant fashion. If we could 
vary ms, we would find that as is reduced from infinity Ui.in1, the critical 
point E in the (T, ^) plane moves toward the /i = axisEJ This is shown 
in Figures 2-4. In nature, E is at some nonzero and ^e- When nig is 
reduced to m^, between Figure 3 and Figure 4, reaches zero. Of course, 



order transition becomes second order. See Figure 1. 




Figure 3: QCD phase diagram for two light quarks and a strange quark with a mass compa- 
rable to that in nature. The presence of the strange quark shifts E to the left, as can be seen 
by comparing with Figure 2. At sufficiently high density, cold quark matter is necessarily in 
the CFL phase in which quarks of all three colors and all three flavors form Cooper pairs. 
The diquark condensate in the CFL phase breaks chiral symmetry, and this phase has the 
same symmetries as baryonic matter which is dense enough that the nucleon and hyperon 
densities are comparable. The unlocking phase transition between the CFL and 2SC phases 
is first order. 



experimentalists cannot vary to^. They can, however, vary fj,. AGS colhsions 
with center of mass priprgy y/g — 5 AGeV create fireballs which freeze out near 



/i - 500 - 600 MeVt2^1l SPS collisions witk^yi = 17 AGeV create fireballs 



which freeze out near n ~ 200 — 300 MeVE23 In time, we will also have data 
from SPS collisions with = 9 AGeV and ^ 12 AGeV and from RHIC 
collisions with y/s = 56, 130 and 200 AGeV and other energies]^ By dialing 
^/s and thus fi, experimenters can find the critical point E. 

We hope that the study of heavy ion collisions will, in the end, lead both 
to a quantitative study of the properties of the quark-giuon plasma phase at 
temperatures well above the transition and to a quantitative understanding of 
how to draw the phase transition region of the phase diagram. Probing the 
partonic matter created early in the collision relies on a suite of signatures 
including: the use of J/^* mesons, charmed mesons, and perhaps the T as 
probes; the energy loss of high momentum partons and consequent effects on 
the high-pT hadron spectrum; and the detection of photons and dileptons over 
and above those emitted in the later hadronic stages of the collision. We will 
not review this program here. Instead, we focus on signatures of the critical 
point. The map of the QCD phase diagram which we have sketched so far is 
simple, coherent and consistent with all we know theoretically; the discovery 
of the critical point would provide an experimental foundation for the central 



"^The firs|r-dntar|from RHIC collisions at y/s = 56 AGcV and y/s = 130 AGeV have already 
appearedE2a'E23 This bodes well for the analyses to come. 
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Figure 4: QCD phase diagram for three quarks which are degenerate in mass and which are 
either massless or Ught. The CFL phase and the baryonic phase have the same symmetries 
and may be continuously connected. The dashed line denotes the critical temperature at 
which baryon-baryon (or quark-quark) pairing vanishes; the region below the dashed line 
is superfluid. Chiral symmetry is broken everywhere below the solid line, which is a first 
order phase transition. The question mark serves to remind us that although no transition 
is required in this region, transition{s) may nevertheless arise as the magnitude of the gap 
increases qualitatively in going from the hypernuclear to the CFL phase. For quark masses 
as in nature, the high density region of the map may be as shown in Figure 3 or may be closer 
to that shown here, albeit with transition(s) in the vicinity of the queptjon mark associated 
with the onset of nonzero hyperon density and the breaking of C/(l)sLj 



qualitative feature of the landscape. This discovery would in addition confirm 
that in higher energy heavy ion collisions and in the big bang, the QCD phase 
transition is a smooth crossover. Furthermore, the discovery of collisions which 
create matter that freezes out near E would imply that conditions above the 
transition existed prior to freezeout, and would thus make it much easier to 
interpret the results of other experiments which study those observables which 
can probe the partonic matter created early in the collision. 

We theorists must clearly do as much as we can to tell experimentalists 
where and how to find E. The "where" question, namely the question of 
predicting the value oi and thus suggesting the ^/s to use to find E, is much 
harder for us to answer. First, as we have stressed in the previous Section, ah 
initio analysis of QCD in its full glory — i.e. lattice calculations — are at 
present impossible at nonzero ^. We must therefore rely on models. Second, 
an intrinsic feature of the picture we have described is that is sensitive 
to the mass of the strange quark, and therefore particularly hard to predict. 
Crude models ^uggest that he could be ~ 600 — 800 MeV in the absence of the 
strange quarkpj'LJ this in turn suggests that in nature ^e may have of order 
half this value, and may therefore be accessible at the SPS if the SPS runs 
with ^/s < 17 AGeV. However, at present theorists cannot predict the value 
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of /is even to within a factor of two. The SPS can search a significant fraction 
of the parameter space; if it does not find E, it wiU then be up to the RHIC 
experiments to map the ^ < 250 MeV region. 

Ahhough we are trying to be helpful with the "where" question, we are 
not very good at answering it quantitatively. This question can only be an- 
swered convincingly by an experimental discovery. What we theorists can do 
reasonably well is to answer the "how" question, thus enabling experimenters 
to answer "where". This is the goal of Ref. 95. The signatures proposed there 
are based on the fact that i? is a genuine thermodynamic singularity at which 
susceptibilities diverge and the order parameter fluctuates on long wavelengths. 
The resulting signatures are nonmonotonic as a function of -y/i: as this control 
parameter is varied, we should see the signatures strengthen and then weaken 
again as the critical point is approached and then passed. 

The critical point E can also be sought by varying control parameters 
other than y^. Ion size, centrality selection and rapidity selection can all be 
varied. The advantage of using ^/s is that we already know (by comparing 
results from the AGS and SPS) that dialing it changes the freeze out chemical 
potential /Lt, which is the goal in a search for E. 

The event-by-event fluctuations of the mean transverse momentum of the 
charged particles in an event, pT, will be enhanced in collisions which freeze 
out near Et^ The fluctuations measured b3j-j!i[A49 at ^/s = 17 AGeV are as 
perfect Gaussians as the data statistics allowli2j as expected for freeze-out from 
a system in thermal equilibrium. The width of the event-by-event distribution^ 
of mean pT is in good agreement with predictions based_an noncritical thermo- 
dynamic fluctuations in an equilibrated resonance gasL3 That is, NA49 data 
are consistent with the hypothesis that almost all the observed event-by-event 
fluctuation in mean p^, an intensive quantity, is thermodynamic in origin. 
This bodes well for the detectability of systematic changes in thermodynamic 
fluctuations near E. 

One analysis described in detail in Ref. 95 is based on the ratio of the 
width of the true event-by-event distribution of the mean pT to the width of the 
distribution in a sample of mixed .events. This ratio was called ^/F. NA49 has 
measured ^fF = 1.002 ± 0.002P^e3 which is consistent with expectations for 
noncritical thermodynamic fluctuations.^ Critical fluctuations of the a field. 



•^This width can be measured even if one observes only two pions per event;|iif| large accep- 
tance data as from NA49 is required in order to learn that the distribution is Gaussian, that 
thermodynamic predictions may be valid, and that the width is therefore the only interesting 
quantity to measure. 

■fin an infinite system made of classical partjefes which is in thermal equilibrium, \/F = 
1. Bose effects increase \/F by 1 — 2%; t^'LJ an anticorrelation introduced by energy 
conservation in a finite system — when one mode fliie*uates up it is more likely for other 
modes to fluct|cra*e down — decreases \/F by 1 — 2%;E3 two-track resolution also decreases 
y/F by 1 — 2%E£3 The contributions due to correlations introduced by resonance decays and 
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i.e. the characteristic long wavelength fluctuations of the order parameter near 
_B, influence pion momenta via the (large) ainr coupling and increase VT'Ej 
The effect is proportional to ■Cfrcczoout' where ^frcczoout isJbhe cr-field correlation 
length of the long-wavelength fluctuations at freezeouto If ^frcozoout 3 fm (a 
reasonable estimate, as we describe below) the ratio ^/F increases by ^.i— 5%, 
ten to twenty times the statistical error in the present measurementlEj This 
observable is valuable because data on it has been analyzed and presented by 
NA49, and it can therefore be used to learn that Pb+Pb collisions at 158 AGeV 
do not freeze out near E. The 3 — 5% nonmonotonic variation in \/F as a 
function of \/s which we predict is easily detectable but is not so large as to 
make one confident of using this alone as a signature of E. 

Once E is located, however, other observables which are more sensitive to 
critical effects will be more useful. For example, a V-Psoft , defined using only 
the softest 10% of the pions in each event, will be much more sensitive to the 
critical long wavekngth fluctuations. The higher px pions are less affected by 
the a fluctuations^ and these relatively unaffected pions dominate the mean 
Pt of all the pions in the event. This is why the increase in ^/F near the critical 
point will be much less than that of V-Fsoft- Depending on the details of the 
cut used to define it, VT^oft should be enhanced by many tens of percent in 
collisions passing near E. Ref. 95 suggests other such observables, and more 
can surely be found. 

The multiplicity of soft pions is an example of an observable which may be 
used to detect the critical fluctuations without an event-by-event analysis. The 
post-freezeout decay of sigma mesons, which are copious and light at freezeout 
near E and which decay subsequently when their mass increases above twice 
the pion mass, should result in a population of pjkins with px ^ 'fn-^jl which 
appears only for freezeout near the critical point lj. If ^froezcout > I/tott, this 
population of unusually low momentum pions will be comparable in number 
to that of the "direct" pions (i.e. those which were pions at freezeout) and 
will result in a large signature. This signature is therefore certainly large for 
^freezeout ~ 3 fm and would not increase much further if ^freezeout were larger 
still. 

The variety of observables which should all vary nonmonotonically with 
i/i (and should all peak at the same -^/s) is sufficiently great that if it were 
to turn out that [le < 200 MeV, making E inaccessible to the SPS, all four 
RHIC experiments could play a role in the study of the critical point. 

The purpose of Ref. 106 is to estimate how large ^freezeout can become, thus 
making the predictions of Ref. 95 for the magnitude of various signatures more 
quantitative. The nonequilibrium dynamics analyzed in Ref. 106 is guaranteed 
to occur in a heavy ion collision which passes near E, even if local thermal 

due to fluctuations in the flow velocity are each much smaUer than 1%0 
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equilibrium is achieved at a higher temperature during the earher evolution 
of the plasma created in the collision. If this plasma were to cool arbitrarily 
slowly, ^ would diverge at Te- However, it would take an infinite time for ^ to 
grow infinitely large. Indeed, near a critical point, the longer the correlation 
length, the longer the equilibration time, and the slower the correlation length 
can grow. This critical slowing down means that the correlation length cannot 
grow sufficiently fast for the system to stay in equilibrium. We use the theory 



of dynamical critical phenomen£E£3 to describe the effects of critical slowing 
down of the long wavelength dynamics near E on the time development of the 
correlation length. The correlation length does not have time to grow as large 
as it would in equilibrium: we find ^froczcout '^/Te ^ 3 fm for trajectories 
passing near E. Although critical slowing down hinders the growth of ^, it also 
slows the decrease of ^ as the system continues to cool below the critical point. 
As a result, ^ does not decrease significantly between the phase transition and 
freezeout. 

We have learned much from the beautiful gaussian event-by-event fluctu- 
ations observed by NA49. The magnitude of these fluctuations are consistent 
with the hypothesis that the hadronic system at freezeout is in approximate 
thermal equilibrium. These and other data show none of the non-gaussian fea- 
tures that would signal that the system had been driven far from equilibrium 
either by a rapid traversal of the transition region or by the bubbling that 
would occur near a strong first order phase transition. There is also no sign 
of the enhanced, but still gaussian, fluctuations which would signal freezeout 
near the critical point E. Combining these observations with the observation 
of tantalizing indications that the matter created in SPS collisions is not well 
described at early times by hadronic model£23 suggests that collisions at the 
SPS may be exploring the crossover region to the left of the critical point E, 
in which the matter is not well-described as a hadron gas but is also not well- 
described as a quark-gluon plasma. This speculation could be confirmed in 
two ways. First, if the SPS is probing the crossover region then the coming ex- 
periments at RHIC may discover direct signatures of an early partonic phase, 
which are well-described by theoretical calculations beginning from an equili- 
brated quark-gluon plasma. Second, if ^/s — 17 AGeV collisions are probing 
the crossover region not far to the left of the critical point E, then SPS data 
taken at lower energies would result in the discovery of E. If, instead, RHIC 
were to discover E with /i_E < 200 MeV, that would indicate that the SPS 
experiments have probed the weakly first order region just to the right of E. 
Regardless, discovering E would take all the speculation out of mapping this 
part of the QCD phase diagram. 

Let us now return to the high density, low temperature regions of the 
phase diagram. In this regime, as we have seen in previous sections, the most 
symmetric starting point with which to begin our analysis is the case of three 
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quarks with degenerate masses. In QCD with — m.u^d as in Figure 4, cold 
dense quark matter is in the CFL phase which has the same symmetries as 
the hyperniiclear matter phase, characterized by a condensate of Cooper pairs 
of baryonsEl Furthermore, many non-universal features of these two phases 
correspondu This raises the possibility that quark matter and baryonic matter 
may be continuously connected^ as shown in Figure 4. 

Both CFL quark matter and hypernuclear matter are superfluids, char- 
acterized by a spontaneously broken global U{1) symmetry associated with 
baryon number. Chiral symmetry is broken in these (or this!) phase. This 
means that as we increase the temperature from zero in Figure 4, there may be 
two distinct phase transitions between the CFL/hypernuclear phase and the 
quark-gluon plasma, as superfluidity is lost and chiral symmetry is restored. 
As hypernuclear matter is heated, we expect superfluidity to be lost at a lower 
temperature than that at which chiral symmetry is restored. At high densi- 
ties, on the other hand, we expect only a single transition because both chiral 
symmetry breaking and superfluidity are caused by the diquark condensate 
which locks color to flavor. Figure 4 is drawn accordingly. The transition at 
which chiral symmetry is restored is first order, for the same reasons as at 
/i = 0: there is no infrared fixed point which could describe a second order 
chiral symmetry restoration transition in QCD with three massless quarksEj 

Nature chooses two light quarks and one middle-weight strange quark, 
rather than three degenerate quarks as in Figure 4. As we have discussed in 
Section 3.2, a nonzero nis weakens those condensates which involve pairing 
between light and strange quarks. The CFL phase requires nonzero (us) and 
{ds) condensates which can only exist if the associated gaps are larger than of 
order m^/2/i. If one paaagines increasing m-s at fixed /z, one finds a first order 
unlocking transition: EjO for larger nis only u and d quarks pair and the 2SC 
phase is obtained. Conversely, as nis is reduced in going from Figure 2 to 3 to 
4, the region occupied-Ui the CFL phase expands to encompass regions with 
smaller and smaller /itj'tlFor any nis ^ oo, the CFL phase is the ground state 
at arbitrarily high densitytj For larger values of nis, there is a 2SC interlude on 
the horizontal axis, in which chiral symmetry is restored, before the CFL phase 
breaks it again at high jdensities. For smaller values of m^, the possibility of 
quark-hadron continuitp as shown in Figure 4 arises. It should be noted that 
when the strange and light quarks are not degenerate, the CFL phase may be 
continuous with a baryonic phase in which the densities of all the nucleons and 
hyperons are comparable; there are, however, phase transition(s) between this 
hypernuclear phase and ordinary nuclear matter made of neutrons and protons 
onlyta 

We can now describe what happens in QCD with m„^d < "is < oo as 
in Figure 3 to CFL quark matter as it is heated. Because the (us) and (ds) 
condensates are smaller than the (ud) condensate, they vanish first. At this 
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temperature, we find a first order unlocking phase transition at whicli chiral 
symmetry is restoredj^ Above this transition, we find the 2SC phase in which 
only red and green u and d quarks pair.[^ The Cooper pairs are ud — du flavor 
singlets and the global flavor symmetry SU{2)l x SU{2)ji is intact. There 
is also an unbroken global symmetry which plays the role of U{1)b- Thus, 
no global symmetries are broken in this 2SC phase. There need therefore be 
no transition between the 2SC and quark-gluon plasma phases in Figure 3 
(or in Figures 1 and 2) because neither phase breaks any global symmetries. 
However, this transition, which is second order in inean field theory, is likely 
first ordes in QCD due to gauge field fluctuationsip at least at high enough 
densitytH; 

Throughout our description of the phase diagram, we have assumed that 
the chemical potentials fi^ = fJ-d = l^-s = fJ- are degenerate. Different values 
of the strange quark rris lead to differences between the Fermi momentum for 
s quarks relative to that for u and d quarks. This analysis neglects electro- 
magnetic and weak effects. In particular, looking ahead to the study of color 
superconductivity in compact stars, it neglects the requirement that the quark 
matter be electrically neutral and in beta equilibrium. If the up, down and 
strange quarks were degenerate, quark matter with /i„ = /i^ = /i^ = would 
be electrically neutral and in beta equilibrium. Instead, because ^ mu,d, 
the quark matter in a compactj^ar must have ^d,s — /^u = A*e > 0, where /Ze is 
the electron chemical potentialE3 In quark matter with average quark chemical 
potential /i ~ 400 — 500 MeV, as may occur in a compact star, reasonable esti- 
mates for /ie are in the few tens of MeV. The presence of a small but nonzero 
fie would certainly have quantitative effects in Figures 1 to 4; we shall see in 
Section 6.6 that it can have qualitative effects also. 



6 Color Superconductivity in Compact Stars 

Our current understanding of the color superconducting state of quark matter 
leads us to believe that it may occur naturally in compact stars. The critical 
temperature Tc below which quark matter is a color superconductor is high 
enough that any quark matter which occurs within neutron stars that are more 
than a few seconds old is in a color superconducting state. In the absence of 
lattice simulations, present theoretical methods are not accurate enough to 
determine whether neutron star cores are made of hadronic matter or quark 

9 That this transition is first order can be seen both via the same argument which domon- 
stratpSjitiiat the unlocking transition between the 2SC and CFL phases at T = is first 
ordefjliZI and via the fact that it is a finite temperature phase transition at which chiral 
symmetry is restored in a three-flavor theory. 

''As we have seen in Section 3, in the 2SC phase at zero temperature, there may be small J = 
1 (ss) and blue (ud) condensates; these condensates would not persist at the temperatures 
we are discussing here. 
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matter. They also cannot determine whether any quark matter which arises 
wiU be in the CFL or 2SC phase: the difference between the u, d and s Fermi 
momenta wih be a few tens of MeV which is comparable to estimates of the gap 
A; the CFL phase occurs when A is large compared to all differences between 
Fermi momenta. Just as the higher temperature regions of the QCD phase 
diagram are being mapped out in heavy ion collisions, we need to learn how 
to use neutron star phenomena to determine whether they feature cores made 
of 2SC quark matter, CFL quark matter or hadronic matter, thus teaching 
us about the high density region of the QCD phase diagram. It is therefore 
important to look for astrophysical consequences of color superconductivity. 



6. 1 Equation of State 

Much of the work on the consequences of quark matter within a compact 
star has focussed on the effects of quark matter on the equation of state, 
and hence on the radius of the star. As a Fermi surface phenomenon, color 
superconductivity has little effect on the equation of state: the pressure is an 
integral over the whole Fermi volume. Color superconductivity modifies the 
equation of state at the ~ (A//x)^ level, typically by a few percentB Such small 
effects can be neglected in present calculations, and for this reason we will not 
attempt to survey the many ways in which observations of neutron stars are 
being used to constrain the equation of stateE 



iir 



We will describe one current idea, however. As a neutron star in a low 
mass X-ray binary (LMXB) is spun up by accretion from its companion, it 
becomes more oblate and its central density decreases. If it contains a quark 
matter core, the volume fraction occupied by this core decreases, the star ex- 



pands, and its moment of inertia increases. This raises the possibilitjtill of a 
period during the spin-up history of an LMXB when the neutron star is gaining 
angular momentum via accretion, but is gaining sufficient moment of inertia 
that its angular frequenqjLis hardly increasing. In their modelling of this effect. 



Glendenning and Webeiiid discover that LMXB's should spend a significant 
fraction of their history with a frequency of around 200 Hz, while their quark 
cores are being spun out of existence, before eventually spinning up to higher 
frequencies. This may exrilain the observation that LMXB frequencies are clus- 
tered around 250-350 HzPi3 which is otherwise puzzling in that it is thought 
that LMXB's provide the link between canonicaLpulsars and millisecond pul- 
sars, which have frequencies as large as 600 HzlH3 It will be interesting to see 
how robust the result of Ref. Ill is to changes in model assumptions and also 
how its predictions fare when compared tp_those of other explanations which 
posit upper bounds on LMXB frequencies tij rather than a most probable fre- 
quency range with no associated upper boundlliil We note here that because 
Glendenning and Weber's effect depends only on the equation of state and not 
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on other properties of quark matter, the fact that the quark matter must in 
fact be a color superconductor wiU not affect the results in any significant way. 
If Glendenning and Weber's explanation for the observed clustering of LMXB 
frequencies proves robust, it would imply that pulsars with lower rotational 
frequencies feature quark matter cores. 

6.2 Cooling by Neutrino Emission 

We turn now to neutron star phenomena which are affected by Fermi surface 
physics. For the first 10^~^ years of its life, the cooling of a neutron star is 
governed by the balance between heat capacity and the loss of heat by neutrino 
emission. How are these quantities affected by the presence of a quark matter 
core? This has been addressed recently in Refs. 115,116, following earlier 
work in Ref. 117. Both the specific heat Cy and the neutrino emission rate 

are dominated by physics within T of the Fermi surface. If, as in the 
CFL phase, all quarks have a gap A ^ T then the contribution of quark 
quasiparticles to Cy and T ,„ is suppressed by ~ exp(— A/T). There may be 
other contributions to L,/lii3 but these are also very pnall. In the CFL phase, 
the specific heat is dominated by that of the electrons c3 (There is an additional 
contribution from the superfluid mode in the CFL phase — i.e. the Goldstone 
boson associated with the spontaneous breaking of U{1)b — and there may 
also be small contributions from the light but not massless pseudo- Goldstone 
bosons associated with chiral symmetry breaking.) Although further work is 
required, it is already clear that both Cy and are much smaller than in the 
nuclear matter outside the quark matter core. This means that the total heat 
capacity and the total neutrino emission rate (and hence the cooling rate) of 
a neutron star with a CFL core will be determined completely by the nuclear 
matter outside the core. The quark matter core is "inert" : with its small heat 
capacity and emission rate it has little influence on the temperature of the star 
as a whole. As the rest of the star emits neutrinos and cools, the core cools 
by conduction, because the electrons keep it in good thermal contact with the 
rest of the star. These qualitative jespectations are nicely borne out in the 
calculations presented by Page et allii3 

The analysis of the cooling history of a neutron star with a quark matter 
core in the 2SC phase is more complicated. The red and green up and down 
quarks pair with a gap many orders of magnitude larger than the temperature, 
which is of order 10 keV, and are therefore inert as described above. Any 
strange quarks present will form an (ss) condensate with angular momentum 
J = 1 -sjjhich locks to color in such a way that rotational invariance is not 
brokencll The resulting gap has been estimated to be of order hundreds of 
keVpZl although applying results of Ref. 49 suggests a somewhat smaller gap, 
around 10 keV. The blue up and down quarks also pair, forming a J = 1 
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condensate which^brcaks rotational invariancei The related gap was estimated 
to be a few keVp but this estimate was not robust and should be revisited 
in light of more recent developments given its importance in the following. 
The critical temperature above which no condensate forms is-|Of order the 
zero-temperature gap A. (T^ = 0.57A for J = condensatesEJ) Therefore, 
if there are quarks for which A ~ T or smaller, these quarks do not pair at 
temperature T. Such quark quasiparticles will radiate neutrinos rapidly (via 
direct URCA reactions like d ^ u + e + i>, u ^ d + + ly, etc.) and the 
quark matter core jdll cool rapidly and determine the cooling history of the 
star as a wholeEiHllij The star will cool rapidly until its interior temperature 
is T < Tc ~ A, at which time the quark matter core will become inert and 
the further cooling history will be dominated by neutrino emission from the 
nuclear matter fraction of the star. If future data were to show that neutron 
stars first cool rapidly (direct URCA) and then cool more slowly, such data 
would allow an estimate of the smallest quark matter gap. We are unlikely to be 
so lucky. The simple observation of rapid cooling would not be an unambiguous 
discovery of quark matter with small gaps; there are other circumstances in 
which the direct URCA processes occur. However, if as data on neutron star 
temperatures improves in coming years the standard cooling scenario proves 
correct, indicating the absence of the direct URCA processes, this would rule 
out the presence of quark matter with gaps in the 10 keV range or smaller. 
The presence of a quark matter core in which all gaps are 3> T can never be 
revealed by an analysis of the cooling history. 



6.3 Supernova Neutrinos 

We now turn from neutrino emission from a neutron star which is many years 
old to that from ttie_protoneutron star during the first seconds of a supernova. 
Carter and ReddjtiS have pointed out that when this protoneutron star is at 
its maximum temperature of order 30-50 MeV, it may have a quark matter 
core which is too hot for color superconductivity. As such a protoneutron star 
core cools over the next few seconds, this quark matter will cool through Tc, 
entering the color superconducting regime of the QCD phase diagram. For 
T ~ Tc, the specific heat rises and the cooling slows. Then, as T drops further 
and A increases to become greater than T, the specific heat drops rapidly. Fur- 
thermore, as the number density of quark quasiparticles becomes suppressed 
by eptp/— A/T), the neutrino transport mean free path rapidly becomes very 
longllis This means that all the neutrinos previously trapped in the now color 
superconducting core are able to escape in a sudden burst. If a terrestrial 
neutrino detector sees thousands of neutrinos from a future supernova, Carter 
and Reddy's results suggest that there may be a signature of the transition 
to color superconductivity present in the time distribution of these neutrinos. 



The Condensed Matter Physics of QCD 



77 



Neutrinos from the core of the protoneutron star will lose energy as they scat- 
ter on their way out, but because they will be the last to reach the surface of 
last scattering, they will be the final neutrinos received at the earth. If they 
are released from the quark matter core in a sudden burst, they may therefore 
result in a bump at late times in the temporal distribution of the detected 
neutrinos. More detailed study remains to be done in order to understand how 
Carter and Reddy's signature, dramatic when the neutrinos escape from the 
core, is processed as the neutrinos traverse the rest of the protoneutron star 
and reach their surface of last scattering. 



6.4 R-mode Instabilities 

Another arena in which color superconductivity comes into play is the physics 
of r-mode instabilities. A neutron star whose angular rotation frequency fl 
is large enough is unstable to the growth of r-mode oscillations which radiate 
away angular momentum via gravitational waves, reducing What does 
"large enough" mean? The answer depends on the damping mechanisms which 
act to prevent the growth of the relevant modes. Both shear viscosity and bulk 
viscosity act to damp the r-modes, preventing them from going unstable. The 
bulk viscosity and the quark contribution to the shear viscosity both becomt 



exponentially small in quark matter with A > T and as a result, as Madseiti^ 
has shown, a compact star made entirely of quark matter with gaps A = 1 MeV 
or greater is unstable if its spin frequency is greater than tens to 100 Hz. Many 
compact stars spin faster than this, and Madsen therefore argues that compact 
stars cannot be strange quark stars unless some quarks remain ungapped. 
Alas, this powerful argument becomes much less powerful in the context of a 
neutron star with a quark matter core. First, the r-mode oscillations have a 
wave form whose amplitude is largest at large radius, outside the core. Second, 
in an ordinary neutron star there is a new source of damping: friction at the 
boundary between the crust and the neutron superfluid "mantle^Ujpe^s the 
r-modes stable regardless of the properties of a quark matter core 
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6.5 Magnetic Field Evolution 

Next, we turn to thcxihysics of magnetic fields within color superconducting 



neutron star cores 00 The interior of a conventional neutron star is a super- 
fluid (because of neutron-neutron pairing) and is an electromagnetic supercon- 
ductor (because of proton-proton pairing) . Ordinary magnetic fields penetrate 
it only in the cores of magnetic flux tubes. A color superconductor behaves 
differently. At first glance, it seems that because a diquark Cooper pair has 
nonzero electric charge, a diquark condensate must exhibit the standard Meiss- 
ner effect, expelling ordinary magnetic fields or restricting them to flux tubes 
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within whose cores the condensate vanishes. This is not the case0 In both the 
2SC and CFL phase, a Unear combination of the U{1) gauge transformation 
of ordinary electromagnetism and one (the eighth) color gauge transformation 
remains unbroken even in the presence of the condensate. This means that 
the ordinary photon and the eighth gluon are replaced by new linear 
combinations 



= cos ao Af, + sin ao , 
A^ = - sin ao + cos ao , (110) 



where A'^ is massless and A-^ is massive. That is, Bq satisfies the ordinary 
Maxwell equations while Bx experiences a Meissner effect, sin(ao) is propor- 
tional to e/j^and turns out to be about 1/20 in the 2SC phase and 1/40 in the 
CFL phasecJ This means that the Q-photon which propagates in color super- 
conducting quark matter is mostly photon with only a small gluon admixture. 
If a color superconducting neutron star core is subjected to an ordinary mag- 
netic field, it will either expel the X component of the flux or restrict it to flux 
tubes, but it can (and doe£2l) admit the great majority of the flux in the form 
of a Bq magnetic field satisfying Maxwell's equations. The decay in time of 

this "free field" (i.e. not in flux tubes) is hmited by the Q-conductivity of the 
quark matter. A color superconductor is not a Q-superconductor — that is the 
whole point — but it may turn out to be a very good Q-conductor due to the 
presence of electrons: if a nonzero density of electrons is required in order to 
maintain charge neutrality, the Bq magnetic field likely decays only on a time 
scale which is much longer than the age of the universec2l This means that a 
quark matter core within a neutron star can serve as an "anchor" for the mag- 
netic field: whereas in ordinary nuclear matter the magnetic flux tubes canJof 
dragged outward by the neutron superfluid vortices as the star spins down 



the magnetic flux within the color superconducting core simply cannot decay. 
Even though this distinction is a qualitative one, it will be difficult to confront 
it with data since what is observed is the total dipole moment of the neutron 
star. A color superconducting core anchors those magnetic flux lines which 
pass through the core, while in a neutron star with no quark matter core the 
entire internal magnetic field can decay over time. In both cases, however, the 
total dipole moment can change since the magnetic flux lines which do not 
pass through the core can move. 



6.6 Crystalline Color Superconductivity and Clitches in Quark Matter 



The flnal consequence of color superconductivity we wish to discuss is the 
possibility that (some) glitches may originate within quark matter regions of 
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a compact star@ In any context in which color superconductivity arises in 
nature, it is hkely to involve pairing between species of quarks with differing 
chemical potentials. If the chemical potential difference is small enough, BCS 
pairing occurs as we have been discussing. If the Fermi surfaces are too far 
apart, no pairing between the species is possible. The transition between the 
BCS and unpaired states as the splitting between Fermi momentaincreases has 



been studied in elm 



HQ superconductorsp3 nuclear superfluid££il and QCD 



superconductorsp30't£j assuming that no other state intervenes. However, 
there is good reason to think that another state can occur. This is the "LQEF" 
state, first explored by Larkin and Ovchinniko'££3 and Fulde and Ferrell££] in 
the context of electron superconductivity in the presence of magnetic impu- 
rities. They found that near the unpairing transition, it is favorable to form 
a state in which the Cooper pairs have nonzero momentum. This is favored 
because it gives rise to a region of phase space where each of the two quarks 
in a pair can be close to its Fermi surface, and such pairs can be created at 
low cost in free energy. Condensates of this sort spontaneously break transla- 
tional and rotational invariance, leading to gaps which vary periodically in a 
crystalline pattern. If in some shell within the quark matter core of a neutron 
star (or within a strange quark star) the quark number densities are such that 
crystalline color superconductivity arises, rotational vortices may be pinned in 
this shell, making it a locus for glitch phenomena. 

The authors of Ref. 49, have explored the range of parameters for which 
crystalline color superconductivity occurs in the QCD phase diagram, upon 
making various simplifying assumptions. We focus primarily on a toy model 
in which the quarks interact via a four-fermion interaction with the quantum 
numbers of single gluon exchange. Also, we only consider pairing between 
u and d quarks, with ij,d — + Sjj, and = /2 — S^i, whereas we expect a 
LOFF state wherever the difference between the Fermi momenta of any two 
quark flavors is near an unpairing transition, including, for example, near the 
unlocking phase transition between the 2SC and CFL phases. 

In the LOFF state, each Cooper pair carries momentum 2q with |q| w 
1.2J/i. The condensate and gap parameter vary in space with wavelength 
7r/|q|. In Ref. 49, we simplify the calculation by assuming that the condensate 
varies in space like a plane wave, leaving the determination of the crystal 
structure of the QCD LOFF phase to future work. We give an ansatz for the 
LOFF wave function, and by variation obtain a gap equation which allows us to 
solve for the gap parameter A^, the free energy and the values of the diquark 
condensates which characterize the LOFF state at a given Sfj, and |q|. We then 
vary |q|, to find the preferred (lowest free energy) LOFF state at a given (5/x, 
and compare the free energy of the LOFF state to that of the BCS state with 
which it competes. The LOFF state is characterized by a gap parameter Aa 
and a diquark condensate, but not by an energy gap in the dispersion relation: 
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we obtain the quasiparticle dispersion relation£j and find that they vary with 
the direction of the momentum, yielding gaps that vary from zero up to a 
maximum of A^. The condensate is dominated by the regions in momentum 
space in which a quark pair with total momentum 2q has both members of the 
pair within ~ A^j of their respective Fermi surfaces. 

The LOFF state is favored for values of 5^ which satisfy < Sfi < dfi2, 
with (5/ii/Ao = 0.707 and 6^2/^0 = 0.754 in the weak coupling limit in which 
Ao <C fi. (For Sfi < Sfii, we have the 2SC phase with gap Aq.) At weak 
coupling, the LOFF gap parameter decreases from 0.23Ao at Sfi — 6^1 (where 
there is a first order BCS-LOFF phase transition) to zero at (5/i = (5/i2 (where 
there is a second order LOFF-normal transition). Except for very close to (5/i2, 
the critical temperature above which the LOFF state melts will be much higher 
than typical neutron star temperatures. At stronger coupling the LOFF gap 
parameter decreases relative to Aq and the window of Sfj,/ Aq within which the 
LOFF state is favored shrinks. The window grows if the interaction is changed 
to weight electric gluon exchange more heavily than magnetic gluon exchange. 

Near the second-order critical point Sfi2, we can describe the phase transi- 
tion with a Ginzburg-Landau effective potential. The order parameter for the 
LOFF-to-normal phase transition is 

<f(r) - -^(ea6ea/33V^""(r)C75^'''(r)) (111) 

so that in the normal phase $(r) = 0, while in the LOFF phase $(r) — 
r^e*^'' "'. (The gap parameter is related to the order parameter by A^ — GTa-) 
Expressing the order parameter in terms of its Fourier modes $(k), we write 
the LOFF free energy (relative to the normal state) as 

Fii^k)}) = J2 (c2(fc2)i$(k)p + aemkt + o(i$n) . (112) 

k 

For (5/i > (5/i2, all of the C2(fc^) are positive and the normal state is stable. 
Just below the critical point, all of the modes $(k) are stable except those 
on the sphere |k| = 2^2, where q2 is the value of |q| at Sfj.2 (so that q2 — 
1.2(5/i2 ~ 0.9Ao at weak coupling). In general, many modes on this sphere 
can become nonzero, giving a condensate with a complex crystal structure. 
We consider the simplest case of a plane wave condensate where only the one 
mode $(k — 2q2) = Fa is nonvanishing. Dropping all other modes, we have 

FiVA) = aiSfx - S^2){Ta)^ + biVA^ , (113) 

where a and b are positive constants. Finding the minimum-energy solution 
for (5/z < (5/X2 , we obtain simple power-law relations for the condensate and the 
free energy: 

rAiSfi)^Kri6fi2~6^i)^/^, F{Sfi)^-KF{6fi2-Sfi)^. (114) 
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These expressions agree well with the numerical results obtained by solving the 
gap equationtJ The Ginzburg-Landau method does not specify the proportion- 
ality factors Ky and Kp, but analytical expressions for these coefficients c^ 
obtained in the weak coupling limit by explicitly solving the gap equation 
yielding 



L28 



GKr = 2V7^Vfe/'5Ai2)2 - 1 ~ I-IS^S;; 



2 (115) 



Notice that because {5fi2^ 6fii)/Sfj,2 is small, the power-law relations ( 114 ) are 
a good model of the system throughout the entire LOFF interval Sfii < Sfi < 
Sfi2 where the LOFF phase is favored over the BCS phase. The Ginzburg- 



Landau expression (113) gives the free energy of the LOFF phase near Sfj,2, 
but it cannot be used to determine the location S/ii of the first-order phase 
transition where the LOFF window terminates. (Locating the first-order point 
requires a comparison of LOFF and BCS free energies.) 

The quark matter which may be present within a compact star will be in 
the crystalhne color superconductor (LOFF) state if 5/i/Ao is in the requisite 
range. For a reasonable value of J/i, say 25 MeV, this occurs if the gap Aq which 
characterizes the uniform color superconductor present at smaller values of (5/i 
is about 40 MeV. This is in the middle of the range of present estimates. Both 
Sfj, and Aq vary as a function of density and hence as a function of radius in 
a compact star. Although it is too early to make quantitative predictions, the 
numbers are such that crystalline color superconducting quark matter may very 
well occur in a range of radii within a compact star. It is therefore worthwhile 
to consider the consequences. 

Many pulsars have been observed to glitch. Glitches are sudden jumps in 
rotation frequency fl which may be as large as A^l/il ~ 10~^, but may also 
be several orders of magnitude smaller. The frequency of observed glitches is 
statistically consistent with the hypothesis that all radio pulsars experience 
glitchestHj Glitches are thought to originate from interactions between the 
rigid neutron star crust, typically somewhat more than a kilometer thick, and 
rotational vortices in a neutron superfluid. The inner kilometer of-crust con- 
sists of a crystal lattice of nuclei immersed in a neutron superfluidE^ Because 
the pulsar is spinning, the neutron superfluid (both within the inner crust and 
deeper inside the star) is threaded with a regular array of rotational vortices. 
As the pulsar's spin gradually slows, these vortices must gradually move out- 
wards since the rotation frequency of a superfluid is proportional to the density 
of vortices. Deep within the star, the vortices are free to move outwards. In the 
crust, however__L[ie vortices are pinned by their interaction with the nuclear 
lattice. Modelfc 
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differ in important respects as to how the stress associated 
with pinned vortices is released in a glitch: for example, the vortices may break 
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and rearrange the crust, or a cluster of vortices may suddenly overcome the 
pinning force and move macroscopically outward, with the sudden decrease in 
the angular momentum of the superfluid within the crust resulting in a sud- 
den increase in angular momentum of the rigid crust itself and hence a glitch. 
All the models agree that the fundamental requirements are the presence of 
rotational vortices in a superfluid and the presence of a rigid structure which 
impedes the motion of vortices and which encompasses enough of the volume 
of the pulsar to contribute significantly to the total moment of inertia. 

Although it is premature to draw quantitative conclusions, it is interest- 
ing to speculate that some glitches may originate deep within a pulsar which 
features a quark matter core, in a region of that core which is in a LOFF 
crystalline color superconductor phase. A three flavor analysis is required to 
estimate over what range of densities LOFF phases may arise, as either {ud), 
{us) or (ds) condensates approach their unpairing transitions. Comparison to 
existing models which descsjibe how p^, p'^ and p^p vary within a quark mat- 
ter core in a neutron staii^ would then permit an estimate of how much the 
LOFF region contributes to the moment of inertia of the pulsar. Furthermore, 
a three flavor analysis is required to determine whether the LOFF phase is a 
superfluid. If thej-anly pairing is between u and d quarks, this 2SC phase is 
not a super|fliiidpl't2l whereas if all three quarks pair in some way, a superfluid 
is obtainedl3ll3 Henceforth, we suppose that the LOFF phase is a superfluid, 
which means that if it occurs within a pulsar it will be threaded by an array 
of rotational vortices. It is reasonable to expect that these vortices will be 
pinned in a LOFF crystal, in which the diquark condensate varies periodically 
in space. Indeed, one of the suggestions for how to look for a LOFF phase 
in terrestrial electron superconductors relies on the fact that the pinning of 
magnetic flux tubes (which, like the rotational vortices of interest to us, have 
normal cores) is expected toJoE much stronger in a LOFF phase than in a 
uniform BCS superconductort23 

A real calculation of the pinning force experienced by a vortex in a crys- 
talline color superconductor must await the determination of the crystal struc- 
ture of the LOFF phase. We can, however, attempt an order o£_magnitude 
estimate along the same lines as that done by Anderson and Itom^ for neu- 
tron vortices in the inner crust of a neulxcux-alar. In that context, this es- 
timate has jSince been made quantitativeE3E23lHi For one specific choice of 
parameterspj the LOFF phase is favored over the normal phase by a free en- 
ergy Floff ^ 5 X (10 MeV)'^ and the spacing between nodes in the LOFF 
crystal is 6 = tt/ (2|q|) ~ 9 fm. The thickness of a rotational vortex is given by 
the correlation length ^ ~ 1/A ~ 25 fm. The pinning energy is the difference 
between the energy of a section of vortex of length b which is centered on a 
node of the LOFF crystal vs. one which is centered on a maximum of the 
LOFF crystal. It is of order Ep ~ Floff ^"^ ~ 4 MeV. The resulting pinning 
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force per unit length of vortex is of order fp ^ E^j}? ^ (4 MeV)/(80 fm^). A 
complete calculation will be challenging because b < and is likely to yield an 
fp which is somewhat less than that we have obtained by dimensional analysis. 
Note that our estimate of fp is quite uncertain both because it is only based 
on dimensional analysis and because the values of A, 5 and -Floff are uncer- 
tain. (We have a good understanding of all the ratios A/Aq, S/j,/ Aq, q/Aq 
and consequently 6Ao in the LOFF phase. It is of course the value of the BCS 
gap Ao which is uncertain.) It is premature to compare our crude result to the 
results of serious calculations of the pinning of crustal neutron vortices as in 
Refs. 135,136,131. It is nevertheless remarkable that they prove to be similar: 
the pinning energy of neutron vortices in the inner crust is i?p « 1 — 3 MeV 
and the pinning force per unit length is /p « (1 - 3 MeV)/ (200 - 400 fm^). 

The reader may be concerned that a glitch deep within the quark matter 
core of a neutron star may not be observable: the vortices within the crystalline 
color superconductor region suddenly unpin and leap outward; this loss of 
angular momentum is compensated by a gain in angular momentum of the 
layer outside the LOFF region; how quickly, then, does this increase in angular 
momentum manifest itself at the surface of the star as a glitch? The important 
point here is that the rotation of any superfluid region within which the vortices 
are able to move_£i:eely is coupled to the rotation of the outer crust on very 



short time scalesEill This rapid coupling, due to electron scattering off vortices 



and the fact that the electron fluid penetrates throughout the star, is usually 
invoked to explain that the core nucleon superfluid speeds up quickly after 
a crustal glitch:__the only long relaxation time is that of the vortices within 



the inner crustliiiil Here, we invoke it to explain that the outer crust speeds 



up rapidly after a LOFF glitch has accelerated the quark matter at the base 
of the nucleon superfluid. After a glitch in the LOFF region, the only long 
relaxation times are those of the vortices in the LOFF region and in the inner 
crust. 

A quantitative theory of glitches originating within quark matter in a 
LOFF phase must await further calculations, in particular a three flavor anal- 
ysis and the determination of the crystal structure of the QCD LOFF phase. 
However, our rough estimate of the pinning force on rotational vortices in a 
LOFF region suggests that this force may be comparable to that on vortices 
in the inner crust of a conventional neutron star. Perhaps, therefore, glitches 
occurring in a region of crystalline color superconducting quark matter may 
yield similar phenomenology to those occurring in the inner crust. This is 
surely strong motivation for further investigation. 

Perhaps the most interesting consequence of these speculations arises in 
the context of xmnpact stars made entirely of strange quark matter. The 



work of WitteiHa and Farhi and Jafff££l raised the possibility that strange 



quark matter may be stable relative to nuclear matter even at zero pressure. 
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If this is the case it raises the question whet 
pulsars, for example — are strange quark stars 



L4C 



£11, observed compact stars — 
rather than neutron stars. 



A conventional neutron star may feature a core made of strange quark matter, 
as we have been discussing above|] Strange quark stars, on the other hand, 
are made (almost) entirely of quark matter with either no hadronic matter 

]t of order one hundred meters thick, which 



content at all or with a thin cm 



contains no neutron superfluidEliHilj The nuclei in this thin crust are supported 



above the quark matter by electrostatic forces; these forces cannot support a 
neutron fluid. Because of the absence of superfluid neutrons, and because of 
the thinness of the crust, no successful models of glitches in the crust of a 
strange quark star have been proposed. Since pulsars are observed to glitch, 
the apparent lack of a glitch mechanism for strange quark star.s has Lean the 
strongest argument that pulsars cannot be strange quark stars M^tNir'^^ This 
conclusion must now bej^yisited. 



Madsen's conclusiortij that a strange quark star is prone to r-mode insta- 
bility due to the absence of damping must also be revisited, since the relevant 
oscillations may be damped within or at the boundary of a crystalline color 
superconductor region. 

The quark matter in a strange quark star, should one exist, would be a 
color superconductor. Depending on the mass of the star, the quark number 



ea^ncrease by a factor of about two to ten in going from the surface to the 
This means that the chemical potential differences among the three 



densit 
center 

quarks will vary also, and there could be a range of radii within which the quark 
matter is in a crystalline color superconductor phase. This raises the possibility 
of glitches in strange quark stars. Because the variation in density with radius 
is gradual, if a shell of LOFF quark matter exists it need not be particularly 
thin. And, we have seen, the pinning forces may be comparable in magnitude 
to those in the inner crust of a conventional neutron star. It has recently been 
suggested (for reasons unrelated to our conaiderations) that certain accreting 



compact starsaaay be strange quark starsliiS although the evidence is far from 
unambiguous £13 In contrast, it has been thought that, because they glitch, 
conventional radio pulsars cannot be strange quark stars. Our work questions 
this assertion by raising the possibility that glitches may originate within a 
layer of quark matter which is in a crystalline color superconducting state. 

There has been much recent progress in our understanding of how the 
presence of color superconducting quark matter in a compact star would affect 
five different phenomena: cooling by neutrino emission, the pattern of the 
arrival times of supernova neutrinos, the evolution of neutron star magnetic 



'Note that a convincing discovery of a quark matter core within an otherwise hadronic 
neutron star would demonstrate conclusively that strange quark matter is not stable at zero 
pressure, thus ruling out the existence of strange quark stars. It is not possible for neutron 
stars with quark matter cores and strange quark stars to both be stable. 



The Condensed Matter Physics of QCD 



85 



fields, r-mode instabilities and glitches. Nevertheless, much theoretical work 
remains to be done before we can make sharp proposals for which astrophysical 
observations can teach us whether compact stars contain quark matter, and if 
so whether it is in the 2SC or CFL phase. 
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